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NOMENCLATURE 

A  Least-squares  constant 
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b.  Lag  window  corresponding  to  no  modification  of  the  power  spectrum 

bp  Harming  lag  window 

B  Least-squares  constant 
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B,  Spectral  window  corresponding  to  b, 

Bp  Spectral  window  corresponding  to  b? 

D  Column  diameter 

D  Particle  diameter 
P 

E  Expectation  value 

f  Frequency 

f„  Nyquist  folding  frequency 

F.  Variance  ratio  with  i  and  j  degrees  of  freedom 

J 

h  Fluidized  bed  height 

h  Packed  bed  height 

k  Equivalent  degrees  of  freedom  of  a  chi -square  distribution 

K  Magnitude  of  the  power  spectrum  of  white  noise  input 

M  Maximum  number  of  time  lags 

N  Total  number  of  data  points 

0  General  output  function 

p  Power  spectrum  of  filtered  data 

P  Modified  power  spectrum 

P,  Power  spectrum  by  B, 
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Pp  Power  spectrum  modified  by  B- 

5  Sum  of  the  squared  errors 

t  Time 

T  Sample  duration 

w.  'actor 

x  General  continuous  random  variable 

x.  Sample  bed  height 

x  Ay     bed  height  (time  average) 

X  Fourier  transform  of  x 

y  Digital  filtered  variable 

Y  Fourier  transform  of  y 

Z  Transfer  function  of  the  system 

Greek  Symbols 

a      :xeak  frequency 

6  Damping  factor 

Y  Prewhitening  filter  coefficient 

r  Transfer  function  of  the  digital  filter 

6  Dirac  delta  function 

n  Integration  varia: 

Autocorrelation  functior 

♦b  Autocorrelation  function  of  band  limited  white  noise 

♦  Power  spectrum 

♦.  Power  spectrum  of  band  limited  white  noise 

♦  Pc     ectrum  of  the  output  of  the  system 
«  Power  spectrum  of  white  noise 
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o  Standard  deviation  of  the  autocorrelation  function 
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t  Correlation  time  constant 
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introduction 

One  of  the  basic  interests  of  the  nuclear  era 
velopment  of  imprcv     ;ctor  types.  In  the  pursuit  of  this  real  I 
vanced  conct :     e  proposed  and  experimentally  tested.  One  such  advanced 
coni     ,ich  has  drawn  considerable  attention  is  the  fluidized  bed  nucL 

jtor  (FBR),  (17).  The  primary  reasons  for  the  Interest  in  this  c 
are     .  ivantares  that  it  has  over  the  fixed  fuel  nuclear  reactors.  Some 
of  the  important  advantages  are  (1*4,17): 

1.  Small  temperature  gradient 

2.  Even  fuel  burn-up  rates 

3.  Low  fuel  fabrication  costs 
k.     Easy  fuel  replacement. 

Previous  experience  with  fixed  fuel  nuclear  reactors  has  shown  that  not  all 
these  properties  are  attainable  in  a  single  reactor.  The  FbR  concept  appears 
to  combine  all  the  properties  in  a  relatively  simple  core  configuration. 
Two  of  the  questionable  properties  of  the  FbR  concept  are: 

1 .  Erosion  and  corrosion  of  the  fuel  elements 

2.  Nuclear  safety  of  a  moving  fuel  nuclear  reactor. 

slon  of  the  fuel  is  due  to  the  friction  caused  when  the  pellets  coll: 
h  one  another  and  with  the  vessel  walls.  The  powder  resulting  is  carried 
in  the  fluid  stream  and  can  cause  damape  as  an  abrasive  and  as  a  result  of 
the  radioactivity  it  carries.  This  problem  has  been  studied  by  Martin  Nuclear 
Company  through  a  contract  with  the  United  States  Atomic  .      Commission. 
Their  experience  with  a  prototype  nuclear  reactor,  the  Oi^anlc  Moder«v 
Fluid! .  .       :tor  (OMFBH),  indicated  that  fuel  erosion  is  a  major  pix>blem 


of*  the  FBR  concept. 

The  nuclear  safety  of  a  system  depends  upon  its  dynamic  behavior.  Nuc- 
lear dynamics  (kinetics)  has  been  studied  extensively  for  a  number  of  year.  . 
As  a  result  one  can  predict  the  behavior  of  a  fixed  fuel  reactor  v/ith  reason- 
able accuracy.  The  problem,  however,  is  compounded  in  a  moving  fuel  reactor 
because  the  core  size  may  be  changing  as  a  function  of  time.  One  must  know 
the  dynamic  behavior  of  the  core  before  he  can  use  the  existing  knowledge  of 
kinetics.  Since  little  is  known  about  time  behavior  of  fluidized  beds,  the 

systems  proposed  to  date  have  been  restricted  to  those  which  have  very 
stable  expansion  behavior  and  a  minimum  of  inhomogeneous  characteristics. 
This  is  a  very  large  restriction  which  might  be  removed  if  more  were  known 
about  the  dynamic  behavior  of  fluidized  beds. 

Various  methods  have  been  devised  to  experimentally  determine  the  dynamic 
behavior  of  the  traditional  "black  box".  They  are  particularly  useful  in  t 
analysis  of  systems  whose  complexity  is  such  that  mathematical  formulation 
is  prohibitive.  The  more  important  of  these  methods  relate  the  output  re- 
sponse to  particular  input  functions,  namely,  sinusoidal,  step,  or  impulse. 
The  sinusoidal  response  method  has  been  applied  successfully  to  fluidiz 

in  the  measurement  of  significant  time  constants  of  the  system  (9). 

Noisi     ,:jis,  another  linportant  method  employed  in  the  :>tudy  c:     :.dc 

I  i.;  used  extensively  in  th     I  of  electro]  I 
Lngful  information  from        in  whJ 

1  to  any  \  i  In  which  1    mtput 

nl  i  I 


:  Lllty  of 

noi.  •  o  riul  No  -uch  application 

The  study  Is  a  preliminary  invest  .  I  .on 

to  determine  w.  or  not  further  work  is  feasible.  A  complt     .  tfjorou>r;i 

study  of  fluldized  bed  dynamics  is  not  intended. 


2.0  TI-IEORETICAL  DEVELOPMENT 

A  complete  development  of  the  statistical  analysis  of  random  functions 
is  beyond  the  scope  of  this  thesis.  Several  good  references  which  describe 
the  application  of  statistical  theory  to  noise  analysis  are  available  (3,5,7, 
8,12).  Chapter  12  of  (5)  is  recommended  as  an  easily  readable  review  of  the 
treatment  of  random  functions,  whereas  (3,7,8,12)  present  a  more  detailed 
analysis.  The  development  which  follows  is  limited  to  the  more  important  con- 
cepts of  noise  analysis. 

The  output  of  the  system  considered  is  assumed  to  be  a  typical  random 
function  from  an  ergodic  stationary  random  process.  The  details  of  this  as- 
sumption are  discussed  in  Appendix  A.  An  important  expression  employed  in 
the  statistical  analysis  of  random  functions  is  the  autocorrelation  function 
defined  as 

*(t)  =  x(t)  x(t  +  t)  (1) 

where  x(t)  is  a  continuous  stationary  random  variable.  x(t  +  t)  is  the  value 
of  the  same  random  variable  sampled  at  a  future  time  determined  by  the  tine 
displacement  variable  x.  The  bar  denotes  an  average.  The  special  case  of 
t  =  0  defines  the  mean  square  value  of  x(t).  If  x(t)  has  a  zero  mean,  <j>(0) 
defines  the  variance  of  x(t).  For  t^  0  the  autocorrelation  function  repre- 
sents the  degree  of  correlation  between  the  value  of  x  at  the  instant  (t)  ana 
the  value  at  some  future  time  (t  +  t). 

The  statistical  average  of  a  random  function  is  generally  defined  as 
nsanble  average.  However,  since  x(t)  is  assumed  to  be  ergodic,  r.q.  (1) 
can  be  expressed  as 


. 


X(t)  x(t  ♦  t)  Ut 

T 

1 


n  of  ' 

n  fund  •  Une  averagi  -  . 

of  the  lnporl     nopertiea  of  <*(0,  the  fact  that      dways  lc:  •  in  or 

equal  to  the  variaj.  ,  .ed  by  cci  Lng  the  inequality 

u  <  |x(t)  1  x(t  +  t)}2  =  x  (t)  +  x2(t  +  t)  t  2x(t)  x(t  ♦  t).   (3) 
Transposir^-,  the  last  term  of  Bq.  (3)  arid  a\       yields 


t   2x(t)  x(t  +  t)  <  x2(l)  +  x2(t  +  t).  (*4) 

Averages  computed  for  stationary  processes     .-.dependent  of  time  dlspl  . 
ment,  thus 


(t  +  t)  =  x2(t)  =  o(0).  (3) 

:  stituting  Eqs.  (5)  and  (1)  into  bq.  CO  yields 

;(t)  <  ?*(0) 
or,  U(t)|  =  ♦(   .  (b) 

The  inequality  in  .  .      •-■comes  ar.     ity  only  when  x(t)  contains  a 
periodic  comjx)nent.   Cf  1  iodic  component  is  2a,  t. 

|*(na)|  =  ♦«)),  t.  = 

The  autocorrelation  function  of  a  to  be  a  series  of 

alt      ly  posit!  U(0)|  spaced  at  into 

of  t  =  na.   I:'  x(t)  is  truly  ran   ,  x(t)  and  x(t  +  t)  become  uncorr- 
for  la         >f  t,  i.e., 


The  autocorrelation  function  can  therefore  be  regarded  as  a  measure  of  t 
randomness  of  a  function. 

Another  useful  property  of  the  autocorrelation  function  is  its  symmetry. 
This  can  be  demonstrated  by  considering 


♦(-t)  = 


lim  1 
T-ko  T 


T 

I   2 


x(t)  x(t  -  t)  dt. 


T 
'2 


Changing  the  variable  to  n  =  t  -  t  yields 


♦(-t)  = 


11m  1 
T-h»  T 


2  "  T 
x(n  +  t)  x(n)  dn. 


"2  "  T 


(7) 


The  interval  of  integration  continues  to  be  2T.  Translation  of  the  time  axis 
does  not  alter  the  statistical  properties  of  a  stationary  process,  therefore 
Eq.  (7)  is  equivalent  to  Eq.  (2),  and 

♦(-t)  =  *(t).  (8) 

It  is  often  convenient  to  express  the  behavior  of  a  function  in  the 
frequency  rather  than  in  the  time  domain.  The  transformation  from  time  to 
frequency  is  typically  performed  by  the  Fourier  transform  defined  as 


X(f)  = 


x(t)  e"iwt  dt 


(9) 


J-, 


where  w  =  2nf 
X(f)  is  called  the  spectral  density  function  of  the  variable  x(t) 

sary  condition  for  the  existence  of  X(f)  is  the  requirement  that 

f  °° 


x(t) 


dt  =  K 


A  neces- 


(10) 


:'orm  is 


x(t)  = 


X(f)  elwt 


nations  (9)  and  (11)  are  commonly  called  r  transform  pal: 

If  x(t)  is  a  stationary  random  function  known  for  all  time,  the  info 

in  Equation  (10)  will  in  general  be  infinite  ai*-  X(f)  will  not      .  Ibei 

fore,  in  order  to  study  the  frequency  cha:  I  sties  of  such  a  function, 


may  consider  a  finite  sepment  of  the  random  function  defined  by 

*p(t) 


m  m 

X(t),   "J  <  t  <  ^ 


(12) 


0, 


t  > 


As  long  as  x(t)  is  bounded,  the  Fourier  transform  of  x^tt)  is  r^iaranteed  to 
exist. 


X^f)      xT(t)  e"iu,t  dt 


(13) 


luivalent  expression  for  the  autocorrelation  function  In  terms  of 
variable  uefined  in  Eq.  (12)  is 


♦  (O  =  T^  T 


T 

(   2 


T 
2 


xT(t)  xT(t  +  t)  dt. 


(1*0 


Due  to  the  definition  of  x^t),  $("0  contains  no  contribution  for  |t|>-,  thus 
Eq.  (1*0  can  be  written  as 

*(t)  =  ^  ^    xT(t)  yt  ♦  t)  dt.  (15) 


Applying  the  Fourier  transform,  Eq.  (9),  to  Eq.  (15)  yields 


♦(f)  = 


f   .  -lux 


{11m  ^ 


T-*» 


T 


x,p(t)  x^t  +  t)  dt}  e  WT  dx 


(16) 


where  <t>(f)  is  the  spectral  density  function  of  4»(t).  If  the  variable  of  inte- 
gration is  changed  to  t  =  n  -  t  and  if  the  order  of  integrating  and  taking  the 
limit  is  interchanged,  Eq.  (16)  becomes 

OO  «   00 


*(f)  =  lim  \ 


Xp(t)  e  w  dt 


/  x  -icon 
x,p(n)  e    dn 


=  lim^(f)  X^f) 


(17) 


T-*» 


where  X^f)  denotes  the  complex  conjugate  of  Xm(f). 

Equation  (17)  is  the  defining  expression  for  the  function  <l>(f),  the  power 
density  spectrum  of  xT,(t).  It  describes  the  square  of  the  frequency  distri- 
bution of  the  random  function  x^/t).  The  name,  power  density  spectrum,  arises 
from  the  original  applications  in  electronics  where  xfp(t)  was  normally  a  cur- 
rent or  voltage.  The  function  x^,(t)  describes  the  power  output  of  a  system 
when  a  resistance  of  one  ohm  is  assumed. 

An  equivalent  expression  can  be  found  for  the  Fourier  transform,  Eq.  (9), 
by  substituting  the  polar  form  of  the  imaginary  exponential, 

-io>t 


=  cos  wt  -  i  sin  cot. 


The  resulting  relation  is 


X(f)  = 


x(t)  cos  wt  dt  -i 


•'-co 


x(t)  sin  oot  dt. 


(18) 


-  *(0  for  x(t)      .  (13)      ♦(f)  acc< 
It  .     "Vlously  shown  that  $(t)  is  an  even  function,  hence  $(t)  sin  wt 
Is  odd  ark:  ;    leal  limits.  The  result 


♦(f)  = 


♦  (t)  cos  on  dr. 


(19) 


It  follows  directly  from  Fourier  analysis  that  the  inverse  of  Eq.  (19)  is 


♦  (t)  = 


♦(f)  cos  wt  df. 


(20) 


The  relation  between  the  power  censity  spectrum  and  the  autocorrelation  func- 
tion exhibited  in  Eqs.  (19)  and  (20)  is  a  statement  of  the  Wiener  theorem  for 
autocorrelation.  Wiener  was  the  first  to  establish  the  validity  of  the  rela- 
tionship for  stationary  random  processes.  A  rigorous  proof  of  the  theoi  • 
riven  by  Lee  (12). 

The  utility  of  the  power  spectrum  of  a  random  function  comes  from  its 
relationship  to  the  transfer  function,  Z(f).  It  is  convenient  to  distinguish 
between  the  amplitude  and  phase  of  Z(f)  by  the  expression 


Z(f)  ■ 


Z(f) 


ie 


(21) 


where  |Z(f)|  is  the  magnitude  and  e  is  the  phase  of  the  transfer  function. 
It  is  shown  in  Appendix  B  that  for  a  given  system 


♦Q(f)  =  %±(t) 


Z(f) 


(22) 


where  ♦  and  ♦.  are  the  power  spectra  of  the  system  output  and  input,  respec- 
tively. The  transfer  function  enters  Eq.  (22)   as  the  square  of  the  magnitude, 
therefore,  according  to  Eq.  (21)  no  phase  information  is  available.  This  is 


10 


a  consequence  of  noise  analysis  based  on  the  autocorrelation  function.  An 
analogous  technique  called  crosscorrelation  may  be  employed  to  determine  the 
transfer  function,  complete  with  phase  information  (12).  The  significant  time 
constants,  however,  can  be  determined  from  the  magnitude  of  the  transfer  func- 
tion. 


Z(f) 


if  both  4>  and  *.  are 
o     1 


Equation  (22)  may  be  used  to  determine 

known.  The  power  spectra  can  be  found  by  measuring  the  input  and  output  of 

the  system  and  applying  Eqs.  (lb)  and  (19).  The  analysis  is  greatly  simplified, 

however,  if  the  input  is  assumed  to  be  white  Gaussian  noise.  Noise  is  said  to 

be  white  if  it  has  a  constant  power  spectrum,  i.e.,  if 

♦  (f)  =  K 
w 

where  K  is  a  constant  and  the  subscript  w  denotes  white  noise.  Rice  (20)  has 
shown  that  white  noise  may  be  regarded  as  a  superposition  of  an  infinite  num- 
ber of  simple  harmonic  oscillators  having  the  following  properties: 

1.  Distribution  continuous  in  frequency 

2.  All  amplitudes  the  same  at  all  frequencies 

3.  Phases  independent  and  random. 

If  the  input  to  a  system  is  white  noise,  then  the  input-output  relation,  Eq. 

(22),  becomes 

2 


♦  (f)  =  K 
o 


Z(f) 


(23) 


Thus,  for'  white  noise  input,  only  the  power  spectrum  of  the  output  is  necessary 


Z(f) 


The  input  to  the  system  may  be  arbitrary  so  long  as  it  i: 


to  determine 
white  noise. 

Klce  has  also  shown  that  a  stationary  Gaussian  random  process  may  be 
regarded  as  the  result  of  passing  whLte  Gaussian  noise  through  a  lim 
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the  system  is  assumed  to  be  linear  and  If  the  output  of  tr*     en 
has  a  normal  distribution  In  steady  state  operation,  the  input  may  be  con- 
sidered as  white  noise. 

It  should  be  noted  that  ideal  white  noise  is  physically  unrealizable  t<  - 
e  its  mean  square  value  is  infinite.  This  can  be  demonstrated  by  consider- 
the  autocorrelation  function  of  white  noise.  Noting  that  the  Fourier-  trans- 
form of  a  constant  is  the  Dirac  delta  function,  6(t),  it  follows  from  Eq.  (20) 


that  $  (t)  is  given  by 


♦u(t)  =  K  6(t). 
w 


Since  6(0)  is  infinite,  the  mean  square  value  of  white  noise,  *w(0)»  Is  also 
seen  to  be  infinite. 

An  approximation  called  band-limited  white  noise  is  physically  realizable. 
This  concept  refers  to  a  random  signal  whose  power  spectrum  is  flat  over  a 
certain  finite  frequency  range  and  zero  elsewhere. 

K,      u1    <    |u)|    <   id~ 

•b(f)    = 

0 ,  |  id  |    <  id,   and   |  id  |    >  iDp 

where  id,  and  id-  are  the  band  limits  of  the  white  noise.  The  mean  square  value 
of  x(t)  is  finite  for  band-limited  white  noise.  This  can  be  exhibited  by  pi 
forming  the  Fourier  transformation  of  ♦.  to  obtain  the  autocorrelation  function, 


♦b(T>- 


=  2 


♦.(f)  cos  idt  df 


r+iD 

K  cos  idt  df 


K  sin  id-t  -  sin  id,  t 

i L  (2*4) 
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The  mean  square  value  is  found  by  evaluating  Eq.  (24)  for  t  =  0.  Applyinj- 
1'  Hospital's  rule  once  gives 

*b<°>  "  I  S  -  ul> 
which  is  finite.  For  practical  applications  noise  may  be  considered  as  white 
if  4>  is  constant  over  the  band  width  of  interest. 


3.0  ANALYSIS  OF  FINITE  AMOUNTS  OF  DATA 

of  the  output  of  a  system  can  be  calculated  from  I 
.  (17)  or  Bq,  (19).  These  two     .  sions  were  shown  to  be  equivalent  in  the 
theoretical  development.  Methods  of  analysis  based  on  both  expressions 

loped  and  applied  successfully  to  the  determination  of  power  spectra 
(1,14,19).  The  more  popular  method  of  estimating  4>(f)  is  based  on  the  Wiener 
theorem,  Lq.  (19).  Since  this  method  depends  on  knowledge  of  the  autocorrela- 
tion function,  it  is  commonly  called  the  autocorrelation  method.  This  method 
is  thoroughly  treated  by  Blackman  and  Tukey  (*0  with  the  applications  to  both 
the  continuous  and  discrete  analysis  of  electrical  sirnals.  Thie  (19)  considers 
the  method  in  its  application  to  nuclear  reactor  noise  analysis.  Both  authors 
treat  the  resolution  and  stability  of  the  calculations  and  develop  relation- 
ships that  enable  one  to  obtain  predetermined  statistical  accuracy  from  experi- 
mental data. 

The  direct  Fourier  analysis  method,  Lq.  (17),  is  treated  by  Akcasu  (1). 
His  expressed  intention  is  to  supplement  the  material  presented  by  Blackman 
and  Tukey  (*0.  He  compared  the  resolution  and  stability  of  the  two  methods  and 
reached  the  conclusion  that  in  general  neither  is  superior  to  the  other.  Al- 
though the  autocorrelation  technique  is  indirect,  the  required  computation 
time  is  no  more  than  the  direct  method.  An  advantage  offered  by  the  autocor- 
relation method  is  the  additional  information  contained  in  the  autocorrelation 
function.  One  quantity  of  importance  which  is  obtained  from  $(t)  Is  the  cor- 
relation time  constant  t  .  This  time  constant  can  be  used  to  estimate  the 

c 

Istical  reliability  of  the  power  spectrum  calculation.  For  this  reason 


in 


the  autocorrelation  method  was  chosen  for  the  analysis  of  the  data  in  this 
thesis.  Details  of  the  method  are  discussed  below. 

The  Wiener  theorem  states  that  the  power  spectrum  of  a  stationary  random 
process  is  related  to  the  autocorrelation  function  of  the  process  by  the  ex- 
pression 


♦  (f)  = 


4>(-r)  cos  o;t  dt 


(19) 


where  <j>(t)  =  lim  ^ 
T-,o> 


T 

t  2 


T 
2 


x^t)  Xp(t  +  t)  dt. 


(15) 


In  the  limit  as  the  sample  duration,  T,  approaches  infinity,  4>(t)  is  defined 
for  all  t.  Under  these  circumstances  *(f)  can  be  calculated  directly  from 

Eq.  (19).  In  practical  applications,  however,  one  deals  with  samples  of  finite 

T 

duration.  The  result  is  a  truncation  of  4>(t)  at  some  maximum,  t  <  -p.  Equa- 
tion (19)  must  therefore  be  modified  to  include  the  effect  of  this  truncation. 
It  is  convenient  to  introduce  a  new  function,  the  lap;  window,  defined  by 


b(t)  <  1  ,  |t |  <  Tm 


bd)  =  0  ,  |t|  >  t 


in 


(25) 


The  product  b(i)  *<{>(t)  c;m   be  used  to  represent  the  truncated  autocorrelation 

function.  Insertion  of  b(-r)  into  Eq.  (19)  modifies  the  expression  as  requ: 
by  truncation  while  retaining  the  identity  of  $(t).  Denoting  the  truncated 
power  spectrum  by  P(f),  Eq.  (16)  becomes 


P(f)  = 


b(i)  4>(t)  e~lu1  ui 


(26) 


The  exponential  form  of  the  Fc 

not  be  symmetrical.  If  d(t)  is      i  as  an  t        Ion, 
)  reduces  to  the  cosine  transform  of  Eq.  (19). 
The  relationship  between  P(f)  and  ♦  (f)  is  found  by  substituting*  the 
Fourier  integral  repn   -  -ion  of  <j»(t)  into  Eq.  (26), 


P(f)  = 


b(T)  e 


-lu>T 


♦  (f)  elu,'T  df  lit. 


J_ 


Rearranging  the  order  of  integration  yields 


P(f)  = 


b(i)  e"i(u)  "  W')T  dx}*(f)  df 


The  quantity  in  brackets  is  recognized  as  the  Fourier  transform  of  b(t).  De- 


fining this  transform  as  D(f)  gives 


P(f)  = 


P(f  -  f  )  ♦(f)  df 


(27) 


■>_» 


Equation  (27)  is  a  convolution  integral  in  the  frequency  domain.  The  w< 

function,  B(f),  is  known  as  the  spectral  window  corresponding  to  the  lag  win- 
dow, b(-r).  The  power  spectrum  computed  from  data  of  finite  duration  is  thus 
seen  to  be  a  weighted  average  of  the  true  power  spectrum.  The  accurate  esti- 
.on  of  ♦(f)  depends  to  a  large  extent  upon  the  frequency  dependence  of  the 
spectral  window. 

pectral  window  is  analogous  to  an  electrical  filter  which  attentat 
all  frequencies  except  those  within  a  certain  band,  Af.  The  more  narrow  the 
band  width  of  a  filter,  the  sharper  the  resolution.  The  same  conditions  apply 
to  the  spectral  window.  The  perfect  spectral  window  is  the  zero  band-width 
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window  described  by  the  Dirac  delta  function. 

B(f  -  f)  =  6(f  -  f) 
Substituting  this  spectral  window  into  Eq.  (27)  yields 

f  °° 


P(f)  = 


6(f  -  f)  *(f)  df'  =  *(f), 


(28) 


i.e.,  perfect  estimation. 

oince  the  delta  function  filter  is  physically  unrealistic,  the  idealized 
band-pass  filter  has  been  introduced.  This  filter  has  unity  transmission 
throughout  the  band  width  and  zero  transmission  elsewhere.  The  frequency 
characteristics  of  the  ideal  band-pass  filter  are  shown  in  Fig.  1.  This 
idealized  model  is  a  reference  by  which  actual  filters  can  be  evaluated. 

The  simplest  lag  window  is  defined  as 


i  .  M  <  t 


m 


b1(x)  = 


0  ,   T 


in 


This  lag  window  corresponds  to  direct  Fourier  transformation  of  the  truncated 
autocorrelation  function.  The  spectral  window  which  b,  (t)  introduces  is  founu 
by  Fourier  transformation. 


B^f)  = 


b,  (t)  cos  u)t  di 


J_ 


ri 


in 


m 


/n  \  sin  2nfi,n 

(1)  cos  cot  dx  =  _jn 

Tlf 


(29) 


The  spectral  estimate  obtained  using  B,(f)  is  found  by  substitutil 
o  :-iq.  (27). 


-4 


o 


o 


• 


0) 


■ 

■ 
o 


• 


ur.\ 


lb 


Px(f)  = 


sin  (2w(f  -  ff)T  ) 

—  *(f')df'.        (30) 

»(f  -  f) 


P,  (f)  is  known  as  the  raw  spectral  estimate  since  it  is  found  from  the  unmodi- 
fled  autocorrelation  function. 

Figure  1  shows  the  frequency  dependence  of  B,(f)  normalized  to  a  band- 
width of  = —  .  It  is  only  a  fair  approximation  to  the  ideal  band-pass  filter. 

m 
The  magnitude  of  the  first  minor  lobe  is  approximately  20%  of  the  major  lobe. 

This  means  that  the  transmitted  power  spectrum  will  contain  significant  con- 
tributions from  frequencies  outside  the  band.  The  first  minor  lobe  contribu- 
tion will  be  negative.  Under  these  conditions  P,  (f)  will  be  a  poor  approxi- 
mation to  *(f),  especially  if  <t>(f)  is  of  significant  magnitude  in  the  vicinity 
of  the  first  minor  lobe. 

A  more  suitable  lag  window,  the  harming  window,  is  defined  by 

|(1  +  COS^),  |l|  <  Tm 

b2(x)  =  m 

0,  |t|   >  T 

*  '  '    m 

The  corresponding  spectral  window  is  found  by  Fourier  transformation  of  b0(t). 


B2(f)  = 


[ 

bp(x)  cos  wt  dr 


P<1  +  cos  — )  cos  uT  dr.  (3D 

'-1  m 

m 


Introducing  the  identity 


cos  A  cos  B  ■  i  cos  (A  +  B)  +  _-  cos  (A  -  B) 
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t     T 


B2(f)  -  ± 


m 


t   T 


COS  u)T  UT  + 


'-i 


I'. 


m 


-T 


cos  2*0'  +  — h  dT 
2t 
m 


m 


♦t 


cos  2n(f  -  ~=-)t 
2t 

m 


(32) 


Substituting';  Lq.  (29)  into  Eq.  (32)  gives 

B2(i')  =  J-  bx(f  -  J-)  +  i  BL(f)  ♦  J-  Mf  ♦  *|-).        (33) 

m  m 

The  spectral  estimate,  Pp(f),  coi -respond inr;  to  the  spectral  window,  bp(f), 

is  thus 


P2(f)  ■  \  P  (f  -  57-)  +  |  P^f)  ♦  \  Px(f  *  27") 

m  m 


(3*0 


where  P,(f)  is  the  raw  spectral  estimate  defined  by  Eq.  (30).  Pp(f)  depends 

on  the  raw  spectral  estimates  spaced  p —  on  either  side  of  the  center  of  the 

m 
banc.   A  logical  choice  for  spacing  between  adjacent  spectral  estimates  is 

also  * — . 
m 

The  frequency  characteristics  of  b~(f)  are  compared  to  b, (f)  in  Fir.  (1). 

The  ramltude  of  the  first  minor  lobe  of  Bp(f)  is  less  than  3.0%  of  the  major 

lobe.  The  transmission  properties  of  bp(f),  tnerefore,  are  better  than  u. 

of  B,(f).  The  resolution  obtained  with  b?(f),  however,  is  less  than  with 

B,(f).  The  equivalent  band-width  of  Bp(f),  i.e.,  the  frequency  band  included 

in  the  major  lobe  of  the  spectral  window,  is  twice  the  band-width  of  b,(f). 

Hence,  spectral  estimates  spaced  by  p —  will  have  considerably  overlapping 

m 
spectral  windows  and  adjacent  estimates  will  not  be  entirely  independent. 
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Alternate  estimates  that  are  separated  by  —  have  less  overlap  of  spectral 

m 
windows,  hence  these  estimates  are  more  nearly  independent. 

Equations  (26)  and  (27)  indicate  two  different  ways  of  employing  the 

lag  window,  b?(i),  to  improve  the  estimation  of  the  power  spectrum.  After 

obtaining  the  truncated  autocorrelation  function  one  may  either 

1.  Modify  $(t)  by  multiplying  bp(t)  and  then 
Fourier  transform  the  product,  or 

2.  Fourier  transform  <t»(x)  to  obtain  P-,(f),  and  then 
convolve  the  product  P-.(f)  B„(f). 

The  two  methods  are  mathematically  equivalent,  however,  the  latter  is  easier 
to  apply.  Equation  (3*0  shows  that  convolution  of  P,  (f)  B„(f)  is  accomplished 

by  smoothing  the  raw  spectral  estimates  spaced  at  intervals  of  a —  with  weights 

111  Tm 

of  ^,  2>  and  }j-. 
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J4.0  DIGITAL  ANALYSIS 

The  theory  developed  thus  far  In  thJ         .lies  directly  to  I 
analysis  of  continuous  samples  of  a  random  function.  If  the  experimental 
data  are  available  only  as  discrete  points  spaced  at  equal  tine  li 
the  analysis  must  be  modified  to  some  extent.  Where  integration  Is  Indicated 
in  the  continuous  case,  summation  or  quadrature  must  be  substituted  in  the 
discrete  case.  These  modifications  are  straightforward  and  will  not  be  dis- 
cussed further  here.  Some  consequences  of  digital  analysis,  however,  are  more 
involved.  These  factors  along  with  methods  of  compensating  for  them  are  dis- 
cussed below. 

An  important  consequence  of  digital  analysis  of  equally  spaced  data 
the  phenomenon  known  as  aliasing.  This  effect  results  from  the  confusion  of 
high  and  low  frequency  content  of  the  sample  function.  The  Nyqulst  sampll: 
theorem  (*0  states  that  a  band-limited  function  can  be  reconstructed  from 
equally  spaced  data  if  the  function  is  sampled  at  two  or  more  points  per  cycle 
of  the  highest  frequency.  As  a  consequence,  any  frequency  higher  than  *-r- 
will  not  be  resolvable  from  data  sampled  at  intervals  of  time,  At.  The  fre- 
quency 

f  -J* 

*N   2  At 

is  known  as  the  Nyquist,  or  folding,  frequency.  High  frequency  information 
is  effectively  folded  back  over  the  low  frequency  spectrum,  I.e.,  all  infor- 
mation resulting  from  frequencies  higher  than  fN  is  detected  within  the  band 
0  <  f  <  f„.  This  result  can  be  shown  by  considering  the  function 

t)  =  sin  2wft. 
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If  g(t)  is  sampled  at  equally  spaced  intervals,  At  =  55-,  it  effectively  be- 
comes 

g(kAt)  =  sin  2irf(kAt),  k  =  0,  1,  2,  ... 
For  multiples  of  the  folding  frequency,  f  =  nf^, 

g(kAt)  =  sin  i2ir(nfM)  (—-) } 

n       ^N 

=  sin  hkv   =  0 
for  all  n  and  k.  The  actual  frequencies  in  the  data  which  are  multiples  of 
fw  cannot  be  distinguished  from  f\,.  Every  frequency  greater  than  f..  will  ap- 
pear within  the  frequency  band  0  <  f  <  f N  under  an  assumed  name,  or  alias. 
The  result  is  an  erroneous  estimate  of  the  power  spectrum.  One  method  of  avoid- 
ing aliasing  effects  is  to  make  certain  that  the  sampling  interval,  At,  is 
chosen  sufficiently  small  so  that  the  folding  frequency  is  higher  than  the 
highest  frequency  of  interest  in  the  power  spectrum. 

A  second  method  of  reducing  the  possibility  of  error  due  to  aliasing  is 
by  filtering  the  data.  Filtering  attenuates  the  frequencies  outside  the  range 
of  interest  and  prevents  their  interference  with  the  analysis.  An  additional 
motive  for  filtering  the  data  is  to  smooth  the  spectrum  so  as  to  make  it  more 
nearly  white.  This  latter  use  of  filtering  is  known  as  prewhitening  and  is 
used  primarily  to  compensate  for  the  non-ideal  nature  of  the  spectral  windows. 
Prewhitening  reduces  the  effect  of  the  contributions  of  the  minor  lobes  by 
flattening  out  the  spectrum.  Filtering  is  typically  carried  out  by  incorpora- 
ting an  electrical  filter  in  the  sampling  circuit.  The  filtering  procedure  can 
also  be  introduced  digitally  after  the  data  has  been  digitalized. 

Digital  filters  are  based  on  a  linear  combination  of  the  digitalized  data. 
A  simple  filter  is  introduced  by  defining  a  new  variable, 
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yJ  =  x.^  +  YXJ-1 


,  is  a  cc.         the  x.  are  the  digital  samples  of  the  or     I  func- 
tion x(t).  Fourier  analysis  of  the  y.  results  in  a  power  spectrum  p(f).  11 
power  spectrum  is  related  to  the  power  spectrum  of  the  oricinal  sample  P(f), 
through  the  transfer  function  of  the  filter. 

In  order  to  obtain  the  relationship  it  is  convenient  to  express  Eq.  (3t>) 
in  an  equivalent  form 

y(t)  =  x(t)  +  Yx(t  -  At). 
Treating  y(t)  as  the  new  input  function  replacing  x(t),  it  is  seen  by  Eq.  (17) 
that  the  new  power  spectrum,  p(f),  is  given  by 


p(f)  = 


lim  1 

T-h»  T 


Y(f) 


Substituting  the  Fourier  integral  representation  for  Y(f)  gives 


p(f)  = 


lim  1 
T-h-  T 


y(t)  e"iwt  dt 


lim  1 


(x(t)  +  Yx(t  - 


At)}  e"lwt  dt 


(36) 


Separating  the  integral  into  a  sum  of  two  integrals  and  substituting 
t  =  t  -  At  in  the  second  term  gives 


/rx  _lim  1 
P(f)  =t—  T 


,.  \   -iut        -iwAt 
x(t)  e    dt  +  ye 


x(x)  e     dt 


lim  1 

T~  T 


X(f)  +  Ye~la)At  X(f) 


P(f) 


1  ♦  ye 


-IwAt 


=  P(f)  r(f). 


(37) 
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Thus,  filtering  by  means  of  a  linear  combination  defined  in  Eq.  (35)  results 
in  multiplying  the  power  spectrum  by  the  transfer  function  of  the  filter, 
r(f).  Expanding  r(f)  gives 


r(f)  = 


2 


1  +  Y  COS  u)At  -  i  y  sin  coAt 

=  1  +  y2  +  2y  cos  wAt.  (38) 

The  desired  function,  P(f),  is  four.d  from  the  expression 

Pm  -  P(f)  -     P(D c>q) 

P{n  ~  TUT'   1  +  y2  +  2Y  cos  oiAt  *  Uyj 

More  complicated  filters  may  be  introduced  by  including  more  terms  in  the 

linear  combination,  Eq.  (35).  The  choice  of  coefficients  and  number  of  terms 

depends  entirely  upon  the  spectrum  being  treated.  In  general,  the  more  peaks 

there  are  in  the  spectrum,  the  more  terms  will  be  necessary  to  prewhiten. 


5.0  I         AL 

Apparatus 

A  photopyaph  of  the  experimental  fluidized  bed  is  shown  In  Kir..  '<-  ■ 
positioning  of  the  movie  came:     the  lir.htinr.  fixtuiv  relative  to  the  t 
section  is  shown  in  FIf;.  3.  A  schematic  diagram  of  the  entire  assembly  is  also 
presented  in  Fir;.  •*. 

A  30-inch-hir.h,  2-inch-diameter  Pyrex  column  was  used  for  the  experimental 
bed  section.  This  section  was  preceded  by  a  24-inch-hlp^,  2-inch-diame: 
calming  section  filled  with  l/'J-lnch-diameter  ceramic  tower  packing  sphen  . 
The  calming  section  provided  a  tortuous  path  for  the  liquid  to  follow  in  order 
to  produce  a  roughly  flat  velocity  profile  at  the  inlet  to  the  bed  section.  A 

ne  mesh  screen  was  used  to  support  the  bed  and  also  to  act  as  a  final  flow 
distributor.  The  screen  was  held  in  place  by  rubber  gaskets  with  the  same 
I.D.  as  the  column  so  as  to  offer  no  restriction  to  the  liquid  flow. 

For  steady  state  operation  the  water  stream  was  pumpea  by  a  1/2  horse- 
power Allis-Chalmers  centrifugal  pump,  Type  SSR.  The  pump  produced  flow  rates 
up  to  10  gallons  per  minute  with  the  by -pass  valve  closed.  The  flow  to  the 
column  was  regulated  by  a  3/^-inch  rJobe  valve  and  was  measured  with  a  Brooks 
rotameter,  size  R-10M-25-3.  A  fifty-five  gallon  drum  was  used  for  the  n 
voir. 

An  auxiliary  step  input  stream  was  also  included  in  the  system.  The 
source  of  water  for  this  section  was  tap  water.  The  flow  rate  was  measured 
with  a  Brooks  rotameter,  size  R-9M-25-2.  .    inputs  were  produced  using  two 
-inch  Gould  solenoid  valves.  One  soienoid  valve,  Type  QR  (normally  open), 
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Figure  2.  Photograph  of  the  experimental  assembly 
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Figure  3.  Positioning  of  the  camera  and  lighting  fixture 
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In  the  by-pass  line.  The  otl.-     ve,  Type  M  (normally  closed),  was  In 
the  line  which  led  to  the  column.  I-oth  valv-  me 

leal  switch  so  that  they  could  be  actuated  simultaneously.  Upon  clo: 
the  circuit,  the  auxiliary  stream  was  Introduced  Into  the  bed  section  as  a 
P  input. 
The  lb  mm.  motion  picture  camera  used  to  obtain  the  data  was  a  Bolex 
H-lb  Rex  with  a  SOM  berthiot  Lytar  lens.  The  focal  length  of  the  lens  was 
25  mm.  and  the  maximum  lens  openinr,  was  f/1.8.  A  Quartz-King  Duo-Light  650 
provided  the  light  source.  A  white  cardboard  background  was  positioned  be- 
hind the  experimental  bed  section  to  give  added  contrast  between  the  bed  and 
the  surroundings.  A  scale  pxaduated  in  tenths  of  Inches  was  included  within 
the  field  of  vision  of  the  camera.  A  Weston  Master  TV  exposure  meter  was  used 
to  measure  the  exposure  settings.  The  Invercone  attachment  was  employed  to 
measure  the  Incident  light.  Eastman  Tri-X  Reversal  Safety  Film,  Type  727b, 
was  used.  The  ASA  rating  of  this  film  is  160. 

5.2  Particles 

The  particles  used  In  the  fluidized  bed  were  copper  particles  obtained 
from  Indiana  Copper  Corporation  (type  S-650).  The  particles  were  obtained  in 
bulk  form  which  contained  various  sizes  and  shapes.  Initial  separation  into 
size  groups  was  accomplished  by  hand  shaking  of  Taylor  Standard  sieves.  A: 
sufficient  particles  were  obtained  in  each  size  group  desired,  the  fina: 
lng  was  carried  out  on  a  shake  table.  The  shaking  operation  was  continued 
for  one  hour  for  each  batch.  The  particles  were  then  graded  for  shape.  T: 
was  accomplished  by  scattering  a  few  particles  on  an  inclined  smooth 
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surface  and  rejecting  all  those  which  did  not  roll  down  the  plane.  This 
procedure  eliminated  the  majority  of  the  ill-shaped  particles  and  left  only 
those  which  were  essentially  spherical.  The  density  of  the  remaining  particles 
was  then  measured  by  a  liquid  displacement  technique.  The  sizes  of  the  par- 
ticles used  in  this  thesis  along  with  their  respective  densities  and  average 
diameters  are  listed  in  Table  I.  The  average  diameters  are  based  on  the  geo- 
metric average  of  the  seive  openings. 

Table  I.  Properties  of  the  bed  particles 


... 

Bed 

Index 

Mesh  Cut 

Density 
(lb/ft3) 

Average  Diameter 
(inches) 

A 

20- 

-24 

530.6 

0.0305 

B 

12- 

■Ik 

527.5 

0.0610 

C 

16- 

-18 

518. 2 

0.0^33 

5 . 3  Procedure 

The  particles  were  loaded  into  the  column  and  a  water  flow  rate  was 
chosen  to  produce  the  desired  bed  expansion  ratio,  h/h  .  Three  expansion 
ratios  were  chosen  for  the  experiment.  These  ratios  were  approximately 
duplicated  for  each  particle  size  in  order  to  keep  h/h  as  a  parameter  of 
the  study.  Camera  settings  were  made  and  motion  pictures  were  taken  of  the 
upper  bed  Interface.  The  camera  was  equipped  with  a  spring  driven  motor 
with  a  capacity  of  approximately  660  frames  per  winding.  During  collection 
of  data,  sub-runs  of  500  frames  each  were  used.  The  number  500  was  cho: 
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to  allow  8  full  .00-foot  i-oll  of  fllr. 

film  speed  would  not  vary  Blgnifli  :«d  of  tl      run. 

A  Recordak  Model  C  microfilm  reader  was  used  to  obtain  I 
from  tl       ped  films.  The  reader  projected  an  image  on  ti.      of  a 

•  glass  screen.  A  st:     edge  was  placed  on  the  face  of  tl 
screen  and  the  data  points  were  read  visually  from  the  scale  1;.       to- 

iph.  The  purpose  of  the  experiment  was  to  follow  one  point  on  the  bed 
intt      as  a  function  of  time.  IXie  to  the  circular  cross  section  of  tl 
experimental  section,  it  was  necessary  to  read  the  data  at  the  side  of  t 
bed  Image  in  the  photograph.  An  attempt  was  made  to  read  the  data  points 
at  the  center  of  the  bed,  however,  it  was  found  that  one  could  not  dif: 
tiate  between  the  motion  of  the  near  and  far  portions  of  the  bed  interface. 
Reading  at  the  side  of  the  bed  avoided  these  undesirable  two  dimensional 
effects.  Some  of  the  runs  were  read  from  both  sides  of  the  bea  in  order  to 
check  the  assumption  that  the  data  were  representative  of  the  entire  bed 
interface.  All  data  points  were  punched  on  IBM  computer  cards  for  analysis 
on  the  liv.-l^lO  digital  computer. 

Step  input  data  were  obtained  for  several  runs  for  which  steady  state 
data  had  already  been  collected.  Step  sizes  were  chosen  to  introduce  small, 
medium,  and  large  perturbations  into  the  inlet  stream  in  order  to  study  the 
effects  of  the  non-linearities  of  the  system.  Readings  were  recorded  first 
for  a  step  increase  in  flow  rate,  then  for  a  step  decrease  in  flow  ra- 
The  flow  rates  of  the  main  and  auxiliary  streams  were  left  the  same  so  that 
the  positive  and  negative  steps  were  of  the  same  size. 
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6.0  DISCUSSION  AND  RESULTS 

It  is  generally  accepted  that  fluidized  beds  have  two  primary  regimes  of 
flow: 

1 )  homogeneous 

2)  slugging  or  bubbling. 

The  first  of  these  is  characterized  by  uniform  bed  density  and  roughly  cons- 
tant interparticle  spacing.  The  upper  fluid-bed  interface  is  well  defined 
and  remains  quite  stationary.  Homogeneous  fluidized  beds  generally  result 
when  both  the  particle  mass  to  fluid  mass  ratio,  p  /Pf,  and  the  particle  dia- 
meter to  bed  diameter  ratio,  D  /D,  are  relatively  low. 

The  slugging  condition  exists  in  most  fluidized  beds  using  gas  as  the 
fluidizing  medium.  In  these  systems  the  mass  ratio  is  quite  high.  Systems 
in  which  the  diameter  ratio  is  high  also  tend  to  operate  in  the  slugging  re- 
gime. The  slugging  fluidized  bed  has  the  appearance  of  a  boiling  fluid  with 
voids  forming  at  the  bottom  of  the  bed  and  rising  to  the  top  of  the  bed  nomin- 
ally intact.  Under  these  conditions  the  upper  fluid-bed  interface  fluctuates 
rapidly . 

Many  fluidized  beds  are  neither  purely  homogeneous  nor  purely  slugging 
but  have  some  characteristics  of  both  types.  The  copper-water  system  studied 
for  this  thesis  is  such  a  system.  The  mass  ratio  of  the  copper-water  system 
is  relatively  high  causing  a  fluctuating  tendency.  The  diameter  ratio,  how- 
ever, was  kept  small  so  that  the  true  slugging  flow  was  not  experienced.  The 
condition  achieved  was  a  relatively  homogeneous  bed  with  a  fluctuating  fluid- 
bed  interface.  Homogeneity  was  desired  in  order  to  remain  within  the  class 
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of  flu     I  beds  r  nucl<     .ctor  appllc        ine  fluctuation 

p  to  a^ly  the  concept  of  nols< 

There  important  parameters  to  be  consider'  of 

a  fluldlzed  bed.  They  a: 

1)  particle  size 

2)  expansion  ratio 

3)  packed  bed  depth 
h)     bed  diameter 

5)  fluldlzing  medium 

6)  particle  density. 

Any  or  all  of  these  parameters  could  have  an  effect  on  the  dynaric  behavior 
of  the  system.  Of  these  six  parameters  only  the  first  three  were  varied  dur- 
ing the  course  of  this  work.  The  latter  three  were  kept  as  constants  of  the 
system.  The  density  of  the  particles  did  vary  somewhat  as  shown  in  Table  I, 
however,  this  variation  should  not  have  affected  the  overall  results.  The 
film  speed  of  the  camera  was  also  varied.  The  film  speed  is  directly  related 
to  the  statistical  analysis  of  the  data  so  discussion  of  this  effect  will  be 
postponed  until  later  in  this  section. 

Table  II  is  a  list  of  the  experimental  steady  state  runs  and  the  values 
of  the  parameters  in  each  case.  Three  particle  sizes  were  used  with  at  least 

three  different  values  of  h/h  for  eacn  particle  size.  The  last  three  runs 

o 

are  replications  of  A-10,  B-l,  and  C-2,   respectively.  The  data  for  these 
runs  were  read  from  the  opposite  side  of  the  bed  in  order  to  determine  if  the 
results  were  representative  of  the  entire  bed  interface.  Runs  C-l,  D-l,  and 
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Table  II.  List  of  runs  and  the  values  of  the  parameters. 


Run  Number 

Particle  Size 
(mesh  cut) 

h 
( in8hes ) 

h/h0 

Camera  Speed 
(fps) 

A-7 

20-24 

3.8 

1.67 

24 

A-8 

20-24 

3.8 

1.68 

32 

A-10 

20-24 

3.8 

1.73 

48 

A-ll 

20-24 

3.8 

1.44 

48 

A-12 

20-24 

3.8 

2.04 

48 

B-l 

12-14 

3.8 

1.78 

48 

B-2 

12-14 

3.8 

2.08 

48 

B-3 

12-14 

3.8 

2.39 

48 

B-4 

12-14 

3.8 

1.47 

48 

B-5 

12-14 

3.8 

1.65 

32 

C-l 

16-18 

3.8 

1.49 

48 

C-2 

16-18 

3.8 

1.73 

48 

C-3 

16-18 

3.8 

1.93 

48 

D-l 

16-18 

5.7 

1.44 

48 

E-l 

16-18 

7.6 

1.44 

48 

A-10-R 

20-24 

3.8 

1.70 

48 

B-l-R 

12-14 

3.8 

1.78 

48 

C-2-R 

16-18 

3.8 

1.74 

48 
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:y   h  for  cc:  ind  cc. 

ctlons  of  the  phot     hie  film  frum  whic: 
Thest     ions  of  film  were  I     !'rom  runs  A-1U,  l-.-l, 
-.  .       rims  were  chosen  for      y  since  their  h/h  valu<       |  proxi- 
mal .:        '.t  was  observed  that  the  magnitude  of  the  l'Juc- 
tuation  of  the  upper  fluici-bed  interface  was  roughly  proportional  to  I 
of  the  particles,  the  smaller  particles  fluctuating  more  than  the  larger  par- 
ticles. The  variance  of  the  experimental  data  which  represents  a  measure  of 
the  fluctuation  was  computed  for  each  run.  The  results  of  these  calculations 
are  included  in  Appenuix  E.  It  was  found  that  the  variance  increased  with 
increasing  h/h  for  a  riven  particle  size,  i.e.,  the  upper  fluid-bed  interface 
became  more  unstable  with  Inert      bed  expansion. 

One  of  the  basic  requirements  for  application  of  noise  analysis  concepts 
to  a  system  is  that  it  be  a  Gaussian  random  process.  The  frequency  distribu- 
tion of  the  experimental  data  for  each  run  was  found  in  order  to  check  the 
validity  of  this  assumption.  The  frequency  distribution  was  obtained  by  count- 

the  number  of  times  a  particular  bed  height  was  observed  in  a  run.  The 
frequency  of  occurrence  was  then  plotted  versus  bed  height.  Figures  6,  7  and 
8  show  representative  curves  obtained  from  the  data.  The  curves  plotted  in 

figures  were  calculated  by  substituting  th€  experimental  mean  and  vari- 
ance into  the  equation  for  the  Gaussian  distribution, 

f(x)  =  — i-  exp  -{(xjlII2} 
f&   o2  2a2 

where  x  is  the  mean  and  o2  is  the  variance.  Vi-e  fit  of  the  experimental  a 

is  that  the  system  approaches  a  Gaussian  random  process. 
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Run  A-10 


Run  B-l 


Run  C-2 


Figure  b.     L-ixcerpts  from  the  photographic  data 
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The  autocorrelation  function  and  pov.er  spectrum  v/ere  calculated  using  the 

lb! 4-14 10  digital  computer  code  described  in  Appendix  C.  Since  the  data  were 

taken  in  digital  form  at  equally  spaced  intervals  of  time,  it  was  necessary 

to  use  digital  approximations  to  obtain  these  functions.  The  autocorrelation 

function,  Eq.  (14),  was  approximated  by  a  normalized  sum  of  lagged  products, 

1  N-n 
N^n  .\       xi  xi+n 
<J>(Tn)  =  IT ,  n  =  0,  1,  2,  ...  M        (40) 

■  Ji  xi  xi 

where  N  is  the  total  number  of  data  points  per  run  and  M  is  the  maximum  num- 
ber of  time  lags,  M  =  x  /At.  Adjacent  data  points  were  used  so  that  Ax  =  At. 
Figures  9,  10,  and  11  show  the  normalized  autocorrelation  function  for  the 
three  frequency  distributions  in  Figs.  6,  7,  and  8,  respectively.  The  results 
of  the  computer  calculations  for  all  the  runs  are  listed  in  Appendix  E. 

The  autocorrelation  function  was  estimated  for  each  sub-run.  The  cor- 
responding sub-run  estimates  were  then  averaged  to  obtain  <j>(t)  for  the  entire 
run.  M  was  chosen  to  give  a  final  power  spectrum  resolution  of  one  radian  per 
second.  The  Nyquist  sampling  theorem  requires  two  samples  per  cycle  for  resolu- 
tion, therefore  M  was  equal  to  one-half  of  the  film  speed,  or  H  =  l/2At. 

The  statistical  accuracy  of  the  autocorrelation  function  has  been  es- 
timated by  Thie  (19).  He  assumed  that  $(i)  could  be  approximated  by 

4>(x)  =  exp  (-  t/tcJ  (41) 

where  x  is  the  correlation  time  constant,  i.e.,  the  time  at  which  the  auto- 
correlation function  is  down  by  a  factor  of  e.  Equation  (41)  is  exact  only 
for  a  Markov  process  (5),  however,  it  is  sufficient  for  purposes  of  error 
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estimation  for  most  Gaussian  random  processes.  The  expression  for  the  stan- 
dard deviation  which  Thie  exhibits  is 

»-  mi 

where  T  is  the  total  time  duration  of  the  sample.  Table  III  lists  the  various 
values  of  T,  t  ,  and  o  for  the  experimental  data. 

Figures  12,  13,  and  14  show  the  power  spectra  corresponding  to  the  auto- 
correlation functions  shown  in  Figs.  9,  10,  and  11.  These  power  spectra  were 
obtained  by  Fourier  transforming  the  autocorrelation  function  as  indicated  by 
Eq.  (19).  The  integration  was  approximated  by  trapezoidal  quadrature  which 

resulted  in  raw  spectral  estimates  spaced  at  o>.  =  s —  .  The  resulting  expres- 

m 
sion  for  the  power  spectrum  was 

M-l 

*(u>n)  =  At  I  $(0)  +  2  I     $(Ti)  cos  un  t1 

♦  ♦<xB)  cos^xj.  (42) 

For  t  =  MAt  and  M  =  ~-r-  the  resolution  becomes 
m  2At 

QA4- 

Aw  =  p-p  =  1.0  radian/sec. 

The  calculated  power  spectra  are  listed  in  Appendix  E. 

Blackman  and  Tukey  (4)  show  that  estimates  of  *  for  a  Gaussian  random 
process  scatter  about  the  true  value  according  to  the  chl-square  distribution. 
For  spectra  that  are  relatively  flat,  i.e.,  those  containing  no  narrow  peaks, 
the  reliability  of  the  measured  values  of  the  power  spectra  depend  only  on 
the  sample  length  and  on  the  number  of  separate  estimates  of  the  spectrum. 
The  reliability  of  each  estimate  may  be  expressed  in  terms  of  k,  the  equiva- 
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Table  III.  Standard  deviation  of  tlie  autocorrelation  function 


Hun  Number  T  c 

(sec)  (sec) 


t  o 

c 


A-7  82.8  0.160  0.044 

A-8  63.1  0.082  0.036 

A-10  41.5  0.106  0.051 

A-ll  41.8  0.090  0.046 

A-12  HI. 6  0.138  0.058 

B-l  41.5  0.101  0.049 

B-2  41.3  0.146  0.059 

B-3  41.3  0.122  0.054 

B-4  iil.O  0.081  0.044 

B-5  63.3  0.131  0.046 

C-l  41.2  0.098  0.049 

C-2  41.3  0.113  0.052 

C-3  41.1  0.147  0.060 

D-l  43.2  0.117  0.052 

B-l  ill. 4  0.141  0.058 

A-10-R  41.4  0.121  0.054 

B-l-R  41.4  0.091  0.047 

C-2-R  41.5  0.110  0.051 


-10- 


V) 

X) 
•H 

•s 


1* 
s 

t — . 
3 


-20- 


w     ( rad/sec ) 
Fip^ire  12.      Power  spectrum  of  Run  A-10 


-10 


5 

■3 


o 
o 


-    - 


u     (rad/sec) 
ire  13.     Power  Bpectrun  of  Hun  b-1 


H3 


to 

r-i 
OJ 
-C 
•H 
O 

o 

-a 


o 

\ 

3 

o 


a»  ( rad/scc ) 
Figure  1*4.   Power  spectrum  of  Run  C-2 


19 


lent  degrees  of  freedom  of  the  ch'.    ire  distribution.  Hie  quantity  k  I 
be  approximated  by  ion 

k  =  (N  -  |)  (|).  (»3) 

Figure  15  gives  the  80%  confidence  limits  for  a  chi-squarv  .     utlon  as 
a  function  of  k.  The  limits  marked  by  the  dotted  lines  show  the  average  sta- 
tistical spread  of  the  experimental  power  spectra  using  a  resolution  of  one 
radian  per  second.  Since  tt>e  error  analysis  was  only  approximate,  the  limits 

0.8*4  «    <  *    <  1.16  ♦ 

ave  =  meas  =      ave 

were  assumed  to  be  representative  of  the  statistical  reliability  for  all  of 

the  runs. 

The  curves  drawn  in  Figs.  12,  13,  and  Ik   are  the  least  squares  fit  of 

a  second  order  transfer  function  to  the  experimental  power  spectrum.  The 

equation  relating  the  power  spectrum  to  the  transfer  function  of  a  system  is 

♦0(f)  =  K  |z(f)|2.  (23) 

A  secDnd  order  transfer  function  is  represented  by  the  expression 

Z(f)  =  - (M) 

a2   +  2aB(Ju))   ♦   (Jw)* 

where  a  is  the  break  frequency  and  6  is  the  damping  factor.  Substituting 

Eq.  (M)  into  Eq.  (23)  gives 

♦  (f)  =  K  1 .  (45) 

a*  +  2a2  u)2(262  -  1)  ♦  mH 

The  experimental  power  spectrum  was  observed  to  tail-off  at  the  high 

frequency  end  of  the  spectrum.  This  behavior  was  also  reported  by  Conn  (6) 

In  the  power  spectrum  of  a  nuclear  reactor,  lie  attributed  the  behavior  to 
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•  component  of  the  spectrum.  This  uncc  I     -d 
component  was  constant  throughout  the  spectrum  so  that  its  ;       could  I 
compensated  for  by  ad    an  arbitrary  constant  to  Bq.  (45) «  Thus,  the 
tion  used  in  fitting  the  data  was 

A 


♦  (f)  = 


• 


a*  +  2a2u,2(26'  -  1)  ♦  *H 
Equation  (16)  has  four  unknown  constants  (a,  6,  A,  and  D)  which  were  to 
be  found  by  fitting  the  equation  to  the  experimental  data.  The  method  of  L 
squares  was  chosen  as  the  criterion  for  the  best  fit.  The  sum  of  the  squared 
errors  is  given  by 

M 

I 

1=1 


{♦   -J  -Bf* 

°i    1 


where 


F  =  a*  +  2a 


<»  l 


(2V 


1)  ♦  M, 


i    -   '  "  "1  ^     "    "i  * 

Differentiating  S  with  respect  to  A,  B,  a,  and  6,  respectively,  gave  the  fol- 
lowing four  normal  equations  to  be  solved 


3A 


=  0  =  2 


^=0  =  2 
3B   U  d 


|^=0=2 
da 


^=0  =  2 
36 


K  -  %    -  B)  <-i  ) 


O 


°i   ri 


A         1aJ  +  W*  (2B2   -  1) 
K  "ft  -Bl  |  A 1 


V 


8a*w  *B 


(17) 
(18) 
(19) 
(50) 


The  first  two  normal  equations,  Eqs.  (17)  and  (10),  are  linear  in  A  and 
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B  and  were  solved  directly  for  A  and  B  in  terms  of  a  and  6.  Weighting  factors 

given  by 

w.  =  l/4> 
i     o± 

were  used  in  the  least  squares  determination  of  A  and  B.  The  weighting  factors 
were  necessary  in  order  to  get  a  reasonable  fit  to  the  high  frequency  end  of 
the  power  spectrum.  Initial  calculations  without  the  weighting  factors  gave 
unsatisfactory  results  with  deviations  of  the  order  of  500$  at  the  high  fre- 
quencies. In  some  instances  the  least  squares  fit  gave  negative  values  at  high 
frequencies.  This  is  unreasonable  since  the  power  spectrum,  which  is  the 
square  of  the  magnitude  of  the  transfer  function,  can  never  be  negative.  In- 
troduction of  the  weighting  factors  caused  the  weighted  error  at  the  high  fre- 
quency end  to  be  of  the  same  magnitude  as  the  error  at  the  low  frequency  end 
of  the  spectrum.  The  result  was  a  reasonable  fit  to  the  experimental  data. 
The  use  of  the  weighting  factors  did  not  significantly  alter  the  values  ob- 
tained for  a  and  B. 

The  last  two  normal  equations,  Eqs.  (^9)  and  (50),  are  non-linear,  there- 
fore they  could  not  be  solved  directly  for  a  and  6.  An  iterative  technique 
based  on  the  Newton-Raphson  method  was  used  to  solve  for  these  two  variables. 
The  details  of  this  method  are  outlined  in  Appendix  D  so  they  will  not  be  dis- 
cussed here.  No  weighting  factors  were  included  in  the  least  squares  deter- 
mination of  a  and  6. 

The  steps  employed  In  the  iterative  method  were: 

1)  Guess  a  and  B  from  the  experimental  spectrum. 

2)  Calculate  A  and  B  from  Eqs.  (H7)  and  (48). 


53 


Calculate  Aa  and  AB  by  the  Newton-Kaphson  method  . 
Eqs.  W)   and  (50). 
^)  Correct  a  and  6. 

5)  Recalculate  A  and  B  and  continue  to  Iterate  until  the 
corrections  were  less  than  sane  predetermined  t  . 
Table  IV  gives  Uie  results  of  these  calculations.  The  computer  output  for 
these  calculations  is  also  listed  in  Appendix  F. 

Only  the  first  portion  of  the  experimental  power  spectrum  was  useu  in 
the  least  squares  fitting  procedure.  The  last  few  points  of  the  power  spec- 
trum were  erratic.  This  behavior  was  attributed  to  inaccuracies  in  reading 
the  data.  Since  the  data  were  obtained  with  an  accuracy  of  ±0.05  inches,  the 
n  frequency  resolution  would  be  expected  to  be  rather  poor.  For  this  rea- 
son the  least  squares  fit  was  carried  out  over  the  first  two-thirds  of  the 
spectrum  only. 

The  results  discussed  up  to  this  point  have  been  for  the  special  case  of 
no  prewhitening  filter  and  a  resolution  of  one  radian  per  second  in  the  power 
spectrum  calculations.  Additional  calculations  were  made  using  a  digital 
filter  and  a  different  resolution  in  order  to  ascertain  the  effects  of  these 
quantities  on  the  final  results  of  the  calculations.  The  results  of  these 
calculations  are  discussed  below. 

All  of  the  runs  were  analyzed  with  a  prewhitening  filter  with  a  value  of 
0.t>  for  the  prewhitening  coefficient.  The  results  of  these  calculations  are 
listed  in  Table  V.  The  use  of  the  filter  resulted  in  quite  different  auto- 
correlation functions,  however,  when  $(t)  was  Fourier  transformed  and  correc- 
ted by  the  transfer  function  of  the  filter,  the  results  were  nearly  the  same 


Table  IV.  Break  frequency  and  damping  factor  of  the  experimental 
data  with  the  corresponding  squared  error 
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Run  Number 

Break  Frequency 

Damping  Factor 

Squared  Error 

A-7 

1.968 

0.6232 

0.0000217 

A-8 

3.235 

0.5184 

0.0003144 

A-10 

2.666 

0.6284 

0.0001161 

A-ll 

3.005 

0.6084 

0.0000792 

A-12 

2.139 

0.6374 

0.0001591 

B-l 

2.928 

0.6925 

0.0002638 

B-2 

2.824 

1.0197 

0.0001631 

B-3 

2.552 

0.7748 

0.0001201 

B-4 

3-353 

0.6466 

0.0000220 

B-5 

2.550 

0.8265 

0.0000294 

C-l 

2.759 

0.5643 

0.0001628 

C-2 

2.491 

0.6101 

0.0001452 

C-3 

2.286 

0.8044 

0.0000626 

D-l 

2.7^5 

0.7818 

0.0001641 

E-l 

2.397 

0.8198 

0.0003907 

A-10-R 

2.552 

0.7370 

0.0000722 

B-l-R 

2.905 

0.5432 

0.0001154 

C-2-R 

2.607 

0.6375 

0.0001611 

v.  Results  of  the  ana 


Hun  Numt 

luency 

Dampinr.  Factor 

Squar     -or 

A-7 

1.877 

0.6430 

0. 000041b 

A-8 

3.112 

0.5286 

0.0012! 

A-10 

2.609 

0.O332 

0.0009379 

A-ll 

2.693 

0.6235 

0.0005258 

A-12 

2.097 

0.6356 

0.0013386 

B-l 

2.866 

0.7176 

0.0015297 

B-2 

2.766 

1.0*133 

0.0011172 

B-3 

2.510 

0.7770 

0.0007068 

B-4 

3.280 

0.6602 

0.0001532 

B-5 

2.159 

0.8340 

0.0001103 

C-l 

2.660 

0.5672 

0.0016677 

C-2 

2.436 

0.6191 

0.0011011 

C-3 

2.228 

0.8027 

0.0005112 

D-l 

2.692 

0.7805 

0.0013054 

&-1 

2.351 

0.8030 

0.0034571 

A-10-. 

2.50*1 

0.7397 

0.0005616 

B-l-R 

2.846 

0.5^96 

0.0007626 

C-2-R 

2.553 

0.6432 

0.0013065 

rt 


as  for  the  non-filtered  case.  All  of  the  runs  taken  at  ^8  frames  per  second 
gave  break  frequencies  and  damping  factors  which  differed  from  the  non-filtered 
values  by  less  than  3%.     These  deviations  were  considered  insignificant  since 
the  statistical  reliability  of  the  power  spectrum  calculations  was  of  the 
order  of  1555.  However,  the  small  deviation  of  a  and  8  from  the  filtered  to 
the  non- filtered  cases  indicated  that  very  little  if  any  aliasing  had  occurred 
in  the  analysis.  The  squared  error  listed  in  Table  V  is  larger  than  that  in 
Table  IV  because  the  magnitude  of  the  power  spectrum  calculated  from  the  pre- 
whitened  data  v/as  greater. 

The  three  runs  taken  at  slower  film  speeds,  A-7,  A-8,  and  b-5,  showed 
more  deviation  when  the  prewhitening  filter  was  applied.  This  indicated  that 
aliasing  could  have  occurred  in  the  analysis  and  that  the  power  spectra  cal- 
culated from  these  runs  were  of  dubious  validity.  Consequently  these  runs  were 
rejected  from  further  consideration. 

The  remaining  runs  were  analyzed  using  a  resolution  of  two  radians  ; 
second.  The  resolution  resulted  in  k  =  266  equivalent  degrees  of  freedom, 
which,  according  to  Fig.  15,  yieided  a  statistical  spread  of 

.(JU  *    <  4>     <  1.10  4 

ave  =  meas  =      ave 

The  statistical  reliability  of  these  calculations  is  somewhat  hotter,  however, 

tht  resolution  was  found  to  be  unsatisfactory  for  further  analysis.  Attem; 

to  utilize  this  data  In  the  least  squares  analysis  to  find  the  break  frequency 

and  damping  factor  resulted  in  an  unstable  condition  which  would  not  converv  , 

Table  VI  lists  the  nine  experimental  runs  for  which  the  parameter  h/h 

was  duplicated  for  each  particle  size.  There  were  three  particle  sizes  and 
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Table  VI:  List  of  n.     -d  in  I        tlcal  analy. 


Particle  Size 


Bed 

ms  ion 

Ratio 

1.47 

1.71 

2.01 

A-ll 

A-10 

A-12 

C-l 

C-2 

C-3 

B-* 

B-l 

B-2 

0.0305 

0.0^33 

0.0610 


three  values  of  bed  expansion  ratio  for  which  this  was  true.  These  runs 

were  statistically  analyzed  to  ascertain  if  Dp  and  h/h  had  a  significant 

effect  on  the  calculated  values  of  a  and  8.  The  statistical  analysis  was 

carried  out  under  the  direction  of  Krause  (11)  and  consisted  of  an  analysis 

of  the  variance  of  the  calculated  data. 

The  first  statistical  analysis  was  the  univariate  analysis  of  variance 

which  treated  a  and  6  as  Independent  variables.  Table  11-3  of  Snedecor  (18) 

is  a  flow  sheet  for  this  analysis  and  outlines  the  method  of  determining  the 

mean  squares  associated  with  the  two-way  table  of  values.  The  variance  ratio, 

F,  which  is  defined  as 

„  _  mean  square  for  treatments 
mean  square  for  discrepance 

was  used  for  the  significance  test. 

Treating  a  as  an  Independent  variable  gave  the  following  values  of  F. 


F2  k   (h/ho)  =  17.5518 


F5  „  (D  )  =  11.0812. 

2,1   p 
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A  standard  table  of  the  F  distribution  shows  that  both  D  and  h/h  were  sig- 

P       ° 

nificant  in  the  observed  values  of  a  at  the  2.5%  level,  i.e.,  there  is  a 
97.5%  probability  that  both  D  and  h/h  had  an  effect  on  the  observed  values 
of  the  break  frequency. 

Treating  6  as  an  independent  variable  gave  the  following  values  of  F. 

F2  jj  (h/ho)  =  4.2161 

Fp  h   (D  )  =  2.3284. 

The  table  of  F  shows  that  D  and  h/h  were  significant  in  the  observed  values 
of  b  at  the  25%  level.  This  indicates  that  there  is  only  a  75%  probability 
that  D  and  h/h  were  responsible  for  the  variation  observed  in  6,  and  that 
there  is  a  25%  probability  that  the  observations  came  from  a  random  sample. 

Since  a  and  B  were  actually  determined  from  the  same  experimental  curve, 
Krause  suggested  the  use  of  a  multivariate  analysis  of  variance.  This  analy- 
sis treated  a  and  6  as  a  coupled  pair  of  observations  and  gave  a  value  of  F 
for  the  pair.  Section  8.9  of  Anderson  (2)  outlines  the  steps  involved  in  the 
multivariate  analysis  of  variance.  The  results  of  this  treatment  were  the 
following  values  of  F. 

F4,6  (h/ho}  =  13'682 


F4,6  <Dp>  =  5'683 


These  results  show  that  D  was  significant  in  the  observed  values  of  a  and  3 
at  the  5%  level,  while  h/h  was  significant  at  the  0.5%  level.  Thus  h/h  was 
still  the  more  important  parameter  over  the  range  covered  by  the  experimental 
work,  however,  the  effect  of  D  was  significant. 
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The  calculated  values  of  a  and  g  using  Lue 

of  0.!  tests.  In  every     the 

significance  level  war-  r.  I  han  for  the  fUtei 

ca;  . 

plot  of  the  calcui  [uencles  as  a  function  of 

bed  expansion  ratio,  h/h  .  The  values  plotted  are  those  corresponding  to 
the  A,  B,  and  C  runs  taken  at  ^8  frames  per  second  and  calculated  without  a 
prewhitening  filter.  The  lines  drawn  in  Kip;.  16  merely  show  the  trends  of 
tie   data  and  are  not  intended  to  be  a  correlation  of  the  data. 

It  is  apparent  that  from  Kip;.  16  the  break  frequency  decreases  as  the 
bed  expansion  ratio  increases.  The  time  constant  of  a  system  is  Inversely 
proportional  to  the  break  frequency,  therefore,  the  system  time  constant  In- 
creases as  h/h  inert    .  .  '.or  was  observed  for  all  three  particle 
sizes.  The  trends  of  the  data  In  PI  .    also  indicated  that  the  break  fre- 
quency changes  more  rapidly  for  small  particles  than  for  large  particles. 

Kigure  17  shows  the  variation  of  the  dam;     ;ctor  as  a  function  of  bed 
expansion  ratio.  Although  *  lite  erratic,  a  general  trend  of  In- 

creasing bed  expansion  ratio  I      rent.  This  would  indicate  that  the  form 
of  the  transfer  function  cL  i  f  bed  expansion.  Since 

8<1.0  represents  an  underdamped  system,  one  woul  d     St  more  overshoot  for 
step  input  responses  when  the  bed  expansion  ratio  is  small  than  when  it  is 
lar, 

Kigure  18  shows  the  effect  that        the  packed  bed  height  has  on  the 
break  frequency  and  damping  factor.  Th     icle  size  and  expansion  ratio 
were  kept  constant  so  that  the  variation  would  be  aue  to  packed  bed  hi 
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The  analysis  of  the  .  dynamic  I     or 

of  the  fluidized  bed  could  be  descrJ  I  of  the  second  on.: 

function,  Kq.  (^).  In  order  to  check  these  results,  a  I      of  step  input 
responses  were  recorded  and  analyzed.  Table  VII  lists  the  conditions  of  the 

input  response  runs.  One  sequence  of  step  response  curves  is  shown  in 
.;.  19  through  2^.  The  positive  step  responses  began  at  conditions  corres- 
ponding to  steady  state  run  B-l  while  the  negative  step  responses  ended  at 
the  B-l  conditions.  The  curves  induced  In  the  figures  are  the  theoretical 
step  responses  calcdated  froni  th    ■  :  size  and  from  the  least  squares  values 
of  a  and  0  for  run  B-l.  Similar  results  were  observed  for  other  particle  sJ 
and  flow  conditio;.  . 

The     etlcal  step  input  response  for  a  second  order  system  with  a 
ictor  .    "  r*an  unit;.  I'erence  (5)  as 

it 
cos  (a  t  -  o)  (51) 

re  0  =  cos"  8 

GUI  I  a  =  asinO. 

r 

The  relationship  between  the  variables  .    .  (51)      idily  seen  in  the 

Argand  diagram,  shown  in  Figure  2*Ja.  In  this  sketch  a   represents  the  magnitude 

of  the  break  frequency  and  a6  represents  the  real  part  of  a.  The  real  part  of 

the  break  frequency  determines  the  rise  time  of  the  response  while  the  lnagin- 


Table  VII.  List  of  step  input  response  runs 


64 


Run  Number 


Corresponding 
Steady  State  Run 


A(h/hQ) 


S-2-0N 
S-2-0ff 

S-3-0n 

S-5-0n 

S-5-0ff 

S-6-0n 

S-6-0ff 

S-7-0n 

S-7-Off 

S-8-0n 

S-8-Cff 

S-9-0n 

S-9-0ff 

S-10-On 

S-10-Off 

S-ll-On 

S-ll-Off 

S-12-On 

S-12-Off 

S-13-On 


A-10 

A-10 

A-10 

B-l 

B-l 

B-l 

B-l 

B-l 

B-l 

13-4 

B-4 

I  -i 

B-l 

C-2 

C-2 

C-2 

C-2 

C-2 

C-2 

C-2 


0.55 
0.5^ 
0.34 
0.38 
0.45 
0.30 
0.30 
0.15 
0.19 
0.23 
0.18 
0.25 
0.27 
0.30 
0.21 
0.10 
0.12 
0.11 
0.08 
0.24 
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Figure  24a.  Relationship  of  variables  In  Eq.  (bl) 


ary  part  determines  the  overshoot.  The  value  of  8  =»  1  represents  the  criti- 
cally damped  condition. 

The  plots  comparing  the  experimental  and  analytical  step  input  responses 
show  that  the  response  can  be  predicted  quite  well  for  small  step  changes 
(for  A(h/h  )  on  the  order  of  0.15).  As  the  step  size  increased,  the  average 
behavior  of  the  bed  interface  deviated  more  and  more  from  the  predicted  be- 
havior. These  results  can  be  explained  as  being  due  to  the  change  of  the  damp- 
ing factor  as  a  function  of  bed  expansion  ration.  As  bed  expansion  increases, 
the  damping  factor  also  increases,  thus  the  overshoot  of  the  step  response 
decreases.  Consequently,  step  responses  calculated  from  the  initial  conditions 
of  the  bed  as  done  here  predict  more  overshoot  than  would  be  expected  experi- 
mentally. 
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It  was  observed  that  the  negative  step  responses  agreed  better,  in 
general,  with  the  analytical  function  than  the  positive  step  responses. 
These  calculations  were  based  on  the  final  state  of  the  bed,  thus  one  would 
expect  the  agreement  to  be  better,  especially  with  regard  to  the  amount  of 
overshoot  in  the  step  response.  Calculations  based  on  a  and  6  for  an  inter- 
mediate value  of  h/h  would  probably  give  better  results.  Calculations  of 
this  type  were  not  attempted  since  the  data  were  too  erratic  to  formulate 
an  adequate  correlation  between  a,  6,  and  h/h  .  The  results  of  the  step 
responses,  however,  were  sufficiently  close  to  indicate  that  the  application 
of  noise  theory  to  the  study  of  the  dynamic  behavior  of  fluidized  beds  is 
valid,  at  least  for  approximate  results. 


.  0  CONCLUSIONS 

The  results  of  these  invest igat ions  nave  led  to  the  conclusic; 
analysis  is  a  valid  method  for  detexmlnJ    e  break  frequency  and  damping  fac- 
tor of  a  copper-water  fluidized  bed  system.  This  conclusion  was  based  on  the 
observations  enumerated  below. 

1.  The  fluctuation  of  the  upper  fluia-bed  interface  had  a  Gaussian  dis- 
tribution about  the  mean,  therefore  it  could  be  considered  as  a 
Gaussian  random  process  to  which  noise  theory  is  applicable. 

2.  Multivariate  analysis  of  the  variance  indicated  that  the  two  experi- 
mental parameters  were  significant  in  their  effect  on  the  final 
results. 

3.  The  copper-water  fluidized  bed  system  was  described  in  terms  of  a 
second  order  transfer  function  with  a  break  frequency  and  damy  I 
factor  which  varied  as  a  function  of  bed  expansion  ratio. 

4.  The  break  frequencies  and  dampinr,  factors  determined  by  least  s 
analysis  of  the  power  spectrum  adequately  described  the  step  input 
response  of  the  system. 

b.  The  break  frequency  of  the  system  decreased  as  a  function  of  bed 
mansion  ratio  for  constant  particle  size. 
The  damping  factor  of  the  system  increased  as  a  function  of  bed 
expansion  ratio  for  constant  particle  size. 

7.  The  break  frequency  decreased  with  packed  bed  height  for  constant 
particle  size  and  constant  bed  expansion  ratio. 

8.  The  damr     .ctor  increased  with  packed  bed  heip.ht  for  cc. 
particle  size  and  constant  bed  expansion  ratio. 
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8.0  SUGGESTIONS  FOR  FURTHER  STUDY 

This  investigation  was  intended  to  establish  the  validity  of  the  applica- 
tion of  noise  theory  to  fluidized  beds.  For  this  purpose  only  a  few  values 
of  bed  expansion  were  chosen  for  each  particle  size.  An  obvious  extension 
of  this  study  would  be  to  include  more  particles  sizes  and  more  bed  expansions 
ratios  along  with  particles  of  different  densities. 

The  results  of  this  study  tended  to  be  erratic  indicating  that  more  data 
were  needed  to  give  more  reliable  results.  The  method  of  data  acquisition 
was  tedious  since  it  depended  upon  visual  reading  of  each  data  point.  Develop- 
ment of  a  method  to  obtain  the  data  either  as  an  electrical  signal  or  as  a 
strip  chart  recording  would  greatly  facilitate  the  analysis.  The  electrical 
signal  would  allow  continuous  analog  computer  analysis  while  the  strip  chart 
recording  would  allow  more  accurate  reading  of  the  digital  data.  Both  methods 
have  been  applied  to  other  systems  with  good  results. 

The  transfer  function  of  a  system  complete  with  phase  information  is 
available  from  cross-correlation  analysis.  This  application  necessitates  a 
random  input  function  v/hose  time  behavior  is  known  so  that  the  input  power 
spectrum  can  be  determined.  An  input  of  randomly  spaced  step  inputs  of  ran- 
dom duration  would  appear  to  offer  some  promise.  This  type  of  input  could 
be  generated  by  using  the  output  of  a  pulse  radiation  detector  counting  a 
radioactive  source  to  trigger  a  solenoid  system  such  as  that  used  in  this 
Investigation. 
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11.0     APPENDICES 


APPENDIX  A 

Stat:     .  Analysis  of  Continuous  Random  Functions 

A  random  function  is  one  whose  future  Is  not  subject  to  precise  \  I 
dictlc:  .    .-.tead  of  beinr,  deterministic,  a  random  function  can  only 
specified  to  the  extent  that  its  probability  density  function,  p(x,t), 
known.  The  probability  density  function  Is  defined  as  that  function  which, 
•i  integrated  between  two  limits,  yields  the  probability  that  the  random 
function,  x(t),  will  be  within  those  two  limits.  Norma lizat ion  of  p(x,t) 
requires  that 

(p(x,t)  dx  =  1. 

The  statistical  average  of  a  random  function,  commonly  denoted  by  x, 
is  the  first  moment  of  the  probability  density  function. 


i-:(x(t)}  = 


x^t)  p(x,t)  dx.  (A-l) 


The  average  so  designated  is  computed  over  an  infinite  ensemble  of  records, 
x.(t),  each  of  finite  length.  This  average  is  also  known  as  the  mean,  :. 
mathematical  expectation,  the  stochastic  average,  and  the  ensemble  a\ 

A  stationary  random  function  is  a  random  function  whose  probability 
density  function  is  time  independent.  The  statistical  average  of  a  station- 
ary random  function  is  thus  independent  of  the  time  at  which  the  ave: 
computed.  Stated  mathematically 


t|x(t  +  T)|  = 


x.(t  ♦  t)  p(x)  dx  =  E  (x(t)} 
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where  t  is  an  arbitrary  time  displacement  variable  and  p(x)  is  the  time 
independent  probability  density  function. 

The  dependence  of  one  continuous  stationary  random  function,  x(t), 
upon  another,  y(t),  is  expressed  in  terms  of  a  Joint  probability  density 
function,  p(x,y).  The  degree  of  correlation  between  the  two  random  func- 
tions is  given  by  the  correlation  function  defined  as 


E{x(t)  y(t)}  = 


x±(t)  yi(t)  p(x,y)  dx  dy.        (A-2) 


The  joint  probability  density  function  is  generally  normalized  so  as  to  make 
E{x(t)  y(t)  }  equal  to  unity  for  perfect  correlation,  zero  for  no  correlation, 
and  minus  one  for  anti-correlation. 

The  autocorrelation  function  of  a  stationary  random  function  is  a  spe- 
cial case  of  Eq.  (A-2).  This  function  describes  the  statistical  correlation 
between  two  values  of  the  same  random  function  displaced  by  an  amount  of  time 
x.  Replacing  y.(t)  by  x,(t  +  t)  gives 


E(x(t)  x(t  +  t)}  = 


x^t)  x^t  +  t)  p(x)  dx  =  R(x).       (A-3) 


Equation  (A-3)  is  the  defining  expression  for  the  ensemble  autocorrelation 
function,  R(t). 

If  the  probability  density  function  is  unknown,  it  becomes  more  con- 
venient to  express  the  statistical  behavior  of  a  continuous  stationary  ran- 
dom function  in  terms  of  its  time  properties.  The  time  average  of  the  func- 
tion x(t)  is  defined  as 


T 
f  2 

x(t)  dt 

T 
2 


where  T  is  the  duration  of  the  sample.  This  limit  may  or  may  not  ex".  I 
depending  upon  I     havlor  of  x(t).  Also,  the  existence  of  x(t)     not 
guarantee  that  it  will  be  equal  to  E(x(t)|.  The  equality  of  time  and  en- 
semble average  holds  only  for  ergodic  random  processes. 

An  ergodic  random  process  is  defined  as  one  for  which  timewise  sampl 
leads  to  the  same  statistical  result  as  ensemble  sampling.  The  ergodic  hy- 
pothesis can  seldom,  if  ever,  be  Justified  formally  (5),  however,  in  most 
physical  problems  a  random  function  may  be  regarded  as  ergodic  if  it  is  time 
stationary.  An  important  consequence  of  the  ergodic  hypothesis  is  that  any 
function  of  an  ergodic  random  function  is  itself  ergodic.  In  particular, 
if  the  random  process  with  sample  function  x(t)  is  ergodic,  then 


xn(t)  =  E{xn(t)} 
I.e.,  the  moments  of  the  probability  density  function  can  be  determined  from 
timewise  sampling  of  the  random  function. 

The  autocorrelation  function  stated  as  a  time  average  is 


,  s   lim  1 
*(t)  "  t—  T 


T 
f  2 

x(t)  x(t  +  t)  dt. 

T 

2 


If  x(t)  is  the  sample  function  of  an  ergodic  random  process,  then 

♦  (t)  =  R(0. 
Thus,  for  an  ergodic  stationary  random  process,  one  can  determine  the  auto- 
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correlation  function  without  knowing  the  probability  density  function  of 
the  random  function  being  correlated. 
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APPENDIX  B 

ut-Out      latlonship  of  the  Power  Spectrum 

The  response,  y(t),  of  a  linear,  constant  parameter  system  to  g 

uts,  x(t),  may  be  described  by  Its  Impulse  response,  h(t),  and 
the  convolution  Integral  (12). 


y(t)  = 


h(i)  x(t  -  t)  dT 


The  function  h(i)  Is  related  to  the  system  transfer  function  by  the  Fourier 
transform.  To  exhibit  this  relationship  consider  the  transfer  function, 
Z(f),  defined  as 

(B-2) 


Z(f)  "  XOT 


where 


Y(f)  = 


y(t)  e-lwt  dt 


and 


X(f)  = 


x(t)  e"iu,t  dt, 


Solvinf;  for  Y(f)  In  Eq.  (B-2)   and  Fourier  transforming  to  obtain  y(t)  gives 


y(t) 


x(f)  e~iu,,t'  df  Z(f)  eiwt  df. 


(B-3) 


The  impulse  response,  h(t),  is  obtained  by  letting  the  Input,  x(t'), 

(B-3)  be  the  Dirac  delta  function  6(f).  Making  these  substitutions 
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h(t)  = 


6(f)  e"iaj,t'  dt»  Z(f)  elu)t  df. 


(B-iJ) 


_oo  '  _oo 


By  the  properties  of  the  delta  function,  the  integral  over  t  is  unity,  thus 
Eq.  (B-4)  reduces  to 


h(t)  = 


Z(f)  eiwt  df 


(B-5) 


which  is  recognized  as  the  Fourier  transform  of  Z(f).  The  inverse  to  Eq. 
(B-5)  is  given  by 


Z(f)  = 


h(t)  e~ia)t  dt. 


(B-6) 


According  to  Eq.  (16)  the  power  spectrum  of  the  output  of  a  linear 
system,  y(t),  may  be  represented  by 


♦0(«  = 


/  v  -iwT 
4>o(t)  e    di 


(B-7) 


J-< 


where 


a  (    \    .  11m  1 


T 
2 

T 

2 


y(t)  y(t  +  t)  dt. 


Substituting  Eq.  (B-l)  for  y(t)  and  y(t  +  t)  in  Eq.  (B-7)  yields 


•0(f)  - 


^-iwr  11m  1 
T-*»  T 


'_. 


h(z)  x(t  -  z)  dz 


T  J_co 

2 


h(n)  x(t  +  t  -  n)  dn  dt  dt. 


(B-8) 
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Rearranging  the  order  of  lnt    I  Ion  In  (B-8)  gi 


Vr) 


■ 


h(z) 


J_c 


h(n)  j^j 


x(t  -  z) 


T 
2 


x(t  +  t  -  n)  dt  dn  dz  di . 


Con:      .,;  only  the  integration  with  respect  to  t,  let  o  =  t  -  z  so  I 


11m  1 
T-~>  T 


T 
I  2 


T 
2 


x(t  -  z)  x(t  +  t  -  n)  dt 


11m  1 
T-h»  T 


x(o)  x(o  +  z  +  t  -  n)  do 


T 
2 


(z  +  t  -  n). 


(b-10) 


Substitution  of  Eq.  (B-10)  into  Kq.  (B-9)  yields 


♦„(»  ■ 


-iuT 


h(z) 


h(n)  *A(z  +  t 


-  n)  dn  dz  dT.      -11) 


^_  i 


Consider  the  integration  with  respect  to  t  in  (B-ll)  and  let  £  ■  z  ♦  t  -  n. 


•0(f)  - 


h(z) 


h(n) 


♦jCO  e"lwU  +  n  "  z)  d^  dn  dz. 


-12) 


'_<*> 


Rearranging  terms  of  Eq.    (B-12)  gives 


V"  ■ 


h(z)  e         dz 


,  ,    »     -icon 
h(n)  e  dn 


^(O  e"1^  dC 


(B-13) 
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By  comparison  with  Eq.  (B-6),  the  integrals  with  respect  to  z  and  £ 
in  Eq.  (B-13)  are  recognized  as  Z*(f)  and  Z(f)  respectively,  where  the  aste- 
risk denotes  complex  conjugation.  The  integral  with  respect  to  £  is  4>,(f). 
Making  these  substitutions  yields 

♦o(f)  =  Z*(f)  Z(f)  *±(f) 


Z(f) 


t^f). 


Thus  it  is  noted  that  the  power  spectrum  of  the  output  is  related  to  the 
power  spectrum  of  the  input  by  the  square  of  the  magnitude  of  the  transfer 
function. 
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APPENDIX  C 

Description  and  Explanation  of  the  IbY-HUO 
Computer  Program  Used  for  Calculation  of  i 
Autocorrelation  Function  and  Power  Spectrum 

This  computer  code  was  written  to  calculate  the  autocorrelation  func- 
tion and  the  power  spectrum  from  the  experimental  data.  The  expressions 
employed  in  the  program  were  discussed  in  the  body  of  this  thesis  and  are 
restated  here  for  convenience.  They  are: 

a)  the  autocorrelation  function 

,  N-n 
N-n  .£,  1  i+n 

N-  I     Xi2 

N  1=1  X 

b)  the  power  spectrum 

M-l 
♦  (u)n)   =  At($(0)  +  2     I     «(tjl)   cos  u»nT1  +  *(tm)   cos  u>nTj,  (^2) 

c)  the  prewhitening  filter 

xi  =  xi+l  "  Y  xi»  (35) 

d)  the  transfer  function  of  the  prewhitening  filter 

r(u>n)  =  1  +  y2  -  2Y  cos  o)p  (38) 

These  equations  were  programmed  in  PDFTTRAN  IV  language.  The  source  program 
is  listed  and  the  logic  diagram  is  shown  in  this  Appendix.  Table  C-l  defines 
the  various  symbols  used  in  the  program. 
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Table  C-l:  Input  data  and  variables  required  for  the 
IBM-l^lO  autocorrelation  function  and 
power  spectrum  computer  program 


Symbol 


Explanation 


A  Raw  spectral  estimates 

AVE  Average  of  input  data 

AZRO  Raw  spectral  estimate  at  zero  frequency 

COR  Correction  factor  for  prewhitening 

DELT  Time  increment  for  input  data 

FREQ  Frequency  at  which  power  spectrum  is  evaluated 

GAMMA  Prewhitening  coefficient 

KLAG  Maximum  number  of  time  lags  to  be  employed 

NPTS  Number  of  data  points  per  sub-run 

NPTST  Total  number  of  data  points  per  run 

NRUN  Number  of  sub-runs  per  run 

NUM  Run  number 

P  Power  spectral  estimates 

PHI  Autocorrelation  coefficients 

PZRO  Power  spectral  estimate  at  zero  frequency 

RET  Reciprocal  time  increment  (film  speed  in  frames  per  second) 

SIGSQ  Variance  for  the  entire  run 

VAR  Variance  for  each  sub-run 

X  Input  data  to  be  correlated 


The  input  to  the  program  was  arranged  as  follows: 
RST  CARD:  FORMAT(2I^>,2F10.iJ,I5) 

1)       Kuii  i     ,  2)  KLAG:  Number  of  power  spectre] 

desired,  3)  RDT:  Film  speed,  U)     GAMMA:  Prewhitening  coefficient, 
and  5)  f  sub-runs  per  run. 

SECOND  CARD:  P0RMAT(I5) 

NPTS:  Number  of  data  points  per  sub-run 
REMAINING  CARDS:  FORMAT (7F1 OM 

X:  Experimental  data  for  a  sub-run  punched  consecutively,  seven  to  a 
card.  (Multiple  sub-runs  were  fed  in  order,  each  sub-run  preceded 
by  the  appropriate  second  card  described  above . ) 
The  output  of  the  program  is  set  up  for  both  card  and  print-out  form. 
The  card  output  is  to  be  used  in  the  least  squares  analysis  program  in  the 
same  order  as  it  comes  from  this  program.  The  first  card  of  the  output  con- 
tains :.   .  IPEST,  RDT,  GAMMA,  SIGSQ,  and  AVE,  respectively.  The  follow! 
cards  (KLAG  of  them)  list  TIME  LAG,  AUTO  FXN,  FREQ,  and  SPECTRUM,  in  that 
order,  punched  one  set  to  a  card.  The  on-line  printer  output  is  similar  to 
the  card  output  with  the  exception  that  the  columns  of  output  have  headi. 
telling  what  they  are. 

Prewhitening  of  the  data  was  accomplished  by  entering  a  value  of  GA.'- 
different  from  zero  in  the  initial  parameter  card.  The  program  automatically 
computes  the  power  spectrum  of  the  prewhitening  filter  and  corrects  the  exper- 
imental power  spectrum.  The  program  is  written  as  a  one  pass  code  and  ace 
only  one  set  of  data  at  a  time.  The  one  run,  however,  may  contain  as  many 
sub-runs  as  desired. 


Pig.  25:  Logic  diagram  for  autocorrelation  function 
and  power  spectrum  computer  program. 


/Start: 

(     Read  in  program] 
V  and  constants 


Zero  the 
mean 


Set: 
ITER=1 
NPTST=0 
PHI=0 


Prewhiten 
data 


no 


Calculate: 
AVE 


y 

Calculate : 

VAH 

Calculate: 

PHI 

r 


Calculate: 

COR 


Calculate: 
P 


Correct: 


Calculate 
Til 


Output : 
TIME,  Pill 
0.  P 


Read: 
X 


End 


. 


THi  •  COMPUTES  THE  AUTOCORRELATION  FUNCTION  AND  PC' 

OF  THI       DATA  WITH  OR  WITHOUT  PRFWHITFNING 
ION  X(6C0)?PHI (50)»J    )  iPl  5<  )  »VAR  (  1C) 

» 2  F 1 0  •  4  •  ] 

.  .  . 

\r\-\    FORMAT (7F1 0.41 

f,FR   NO  DATA  PTS  F  I  l  M  CPFFD   GAMMA»5X tSHVAR I Al 

1»7X»4H" 

■  '  (5XtI2*9Xtl4t8X»F4«0t5X»F4«2t5X»El0*4t5XtFl0«71 

(//l  KtSHTlME  LAGt7X.8HAUTC  FXN  ♦  1  OX  •QHFREOUEnC  Y  ,  7X  .  3HSPECTRU 
1MJ 

204  F     '  '  (  1  )XtF8.5.5X.Fl  l.f,t  10X,F8.4t5X.FU.o) 

ITFP=1 

NPTST=^ 

PFAP(  1  ,  \n\   )NUM,KLAG,RDT  » GAMMA ,NRUN 
1  1=1.50 
1  PHI ( I )=r. 
?  RFAnil,l^?)NPTc 

RFAD(1.103)(X( II»I«1 »NPTS) 

IF (GAMMA. FO«0«0)G0    TO   * 

DC    3     I=?,NPTS 

J=I-1 

3  X ( J)=X( I )-GAMMA#X( J) 
NPTS=NPTS-1 

4  sir 

DC    5     1=1 .NPTS 

5  SUM-SUM+XJ I ) 
AVE*SUM/FLCAT(NPTS) 

:  f  [*1»NPTS 

6  x  (1  )  =  X{ I )-avc 
Si 

DC  7  1=1 ,NPT« 

7  SUM-SUM+X ( I )*X( I ) 

VAR (  ITFR) =SUM/FLCAT(NPTS) 

DC  9  1  =  1  ,KLAG 

K        I 

Slr 

DC  8  J=1  »| 

L=I  +J 

8  5UI     "  •  *  (  J)  *X  (L  ) 

O  PHI  (  I  )=PHI  (  I  I+«:ijv/IV»RI  ITFR)*FLCAT(NRUN*K  )  ) 
r<PTST  =  NPTST  +  NP' 

I F( •   N.FO.ITER »GC  TC  10 

I  T-Dr  ITFD+1 

GC  T~  ? 
10  S!GSO«O,0 

I  =1 ,NRUN 
20  SIGSO=  :     fV-R(  I  )/FLCAT(NRUN) 

PI « 3. 14159265 
?1  ISUM»KLAG-1 

DFLT=1 ./RDT 

SI  '-  . 
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11 


1? 

13 


14 


15 

16 


17 


DO 

SUM 

AZR 

DC 

SUM 

DC 

ARG 

SUM 

A(  I 

PZR 

P(  1 

L=K 

P(< 

DC 

P(  I 

IF( 

TFP 

PZP 

P  (K 

DC 

CCR 

P(  I 

WRI 

WRI 

WRI 

WRI 

DJM 

DUM 

WRI 

WRI 

do 

TIM 
FRE 
WRI 
WRI 
S  TC 
END 


11  1  =  1 
=SUM+2 
C=DELT 

13  1  =  1 

12  J=l 
=FLCAT 
=SUM+C 
)=DFLT 
0  =  0.5* 
)=0.25 
LAG-1 
LAG)=0 

14  1=2 
)=0.25 
GAMMA. 
M=l .+G 
C=PZRC 
LAG)=P 

15  1  =  1 
=TERM- 
)=P( I ) 
TE(3,2 
TE(2.2 
TE(3,2 
TF(3.2 
=  0.0 
M=]  .0 
TF(3,2 
T^(2.2 
17  1  =  ] 
F=FLOA 
0=FLCA 
TE(3.2 
T  E  ♦  2  ,  2 
P 


>  I  SUM 

.*PHI ( I ) 

*( 1 .+SUM+PHI (KLAG ) ) 

,KLAG 

»  I  SUM 

(  I*J)/F|_OAT(KLAG) 

CS( ARG*PI )*PHI ( J) *2. 

♦d.+SUM  +  CCS  (FLOAT  (  I  )  *P  I  )  *PH  I  (  KL  AG  )  ) 

AZRC+0.5*A( 1 ) 

*AZRC+0.5*A(1)+0.25*A(2) 


.5* 

♦  L 
*A( 
FQ. 
A  MM 
/(T 
(KL 

♦  IS 
2.* 
/CO 
01) 
02) 
02) 
03) 


A(L)+0.5*A(KLAG) 

1-1  )+0.5*A(  I  )+0.25*A(  1+1 ) 

0.0)G0  TO  16 

A*GAMMA 

ERM-2«*GAMMA) 

AG) /(TFRM+2.*GAMMA ) 

UM 

GAMMA*CCS( FLOAT ( I  )*PI /FLOAT (KLAG) ) 

R 

NUM.NPTST,RDT,GAMMA»SIGSQ,AVE 
NUM.NPTST.RDT. GAMMA »S IG5G.AVE 


04 )DUM,DUMM  »DUM,PZPC 

04  )  DUM,PUMM,DlJM,PZPC 

»KLAG 

T( I )/RDT 

T(  I  )  *RDT/FLCAT(2*KLAG) 

04) TIME. PHI  (  I  )  .FRFQ,P(  I  ) 

04) TIME. PHI ( I ) ♦FRFO,P( I ) 


APPENDIX  D 

Description  and     ination  of  the  IbM-HQO 
Conputer  Program  Used  for  Calculation  of  the 
Break  Frequency  and  Damping  Factor 

This  computer  code  was  written  to  calculate  the  break  frequency  and 
damping  factor  by  simultaneous  solution  of  Eqs.  (^7)  through  (50).  Bqua- 
tlons  (*47)  and  (48)  were  solved  for  A  and  B  in  terms  of  a  and  6.  Equations 
(49)  and  (50)  were  then  solved  for  a  and  8  by  the  Newton-Raphson  iteration 
method.  This  code  was  written  to  accept  the  output  of  the  program  described 
in  Appendix  C  directly. 

Equations  (49)  and  (50)  were  denoted  by 

?|-  R^o.B)  =  0  CM) 

§|-  R2(a,B)  =  0.  (D-2) 

Although  A  and  B  were  unknown  quantities,  they  were  considered  as  known 
values  in  the  solution  of  Eqs.  (D-l)  and  (D-2).  Their  values  were  calculated 
from  the  current  values  of  a  and  8  using  Eqs.  (47)  and  (48). 

Expanding  R,  and  Ft,  in  a  Taylor  series  expansion  about  the  initial  guess 
(a  ,6  )  gives 

R^o.6)  -  0  *  R1(.0.8o)  ♦  3-^y°0.B0)  Aa  ♦  ^0Sn)   46     (D-3) 

R2(a  ,6)  -  0  *  R2(a0>60)  ♦  ,-&,<.„  .6  „>  A.  ♦  ^^.a,,)  48     <D-*> 

Solving  Eqs.   (D-3)  and  (D-4)   for    Aaand    AB yields 


9^ 


where 


Aa  -  ^ 


36 

-R   ^2 
36 


D  = 


A6  = 


D 


3R 


1  -R, 


9a 
3R 


2  -Ft 


8a 


9R1  aRj^ 

3a"~  36~ 

3R2  8R2 

~3a  "3T~ 

it  being  understood  that  the  elements  of  the  determinants  are  to  be  evaluated 

at  (a  B  ). 

°  3R.     3R. 

The  derivatives,  -=—  and  -rr.  were  determined  using  the  following  numeri- 

da        a  p 


cal  approximations: 


*R±   ^  R1(1.05ao,3Q)  -*±(*0,»0) 
3a  '  -05a 


(D-6) 


3R 


i  .  Vv1-05^  -  WV 


36 


(D-7) 


.056, 


The  terms  of  the  approximate  derivatives  were  denoted  by  the  matrix  elements 

GlfRl(VV     G12  =  R2(ao'6o) 

G21  =  Rl(ao'1«056o)  G22  =  R2(ao>1'056o) 
G31  =  R1(1.05ao,6o)  G32  =  R2(1.05aQ, 6Q) . 

Making  these  substitutions,  Eqs.  (D-6)  and  (D-7)  become 

Q  *   ^2  x  G21  "  Gll     „   ^1  x  G31  ~  Gll 
21   '   3a  "    .05a    °31     3a  '     .05a 


3R?   G?2  -  G,„ 
22     36      .056     32     36 


^1  .  G32  '  G12 
.056^ 


(D-8) 


where  the  asterisk  ii  ote  the  redefinition  of  the  G  . 

:>ol  was  used  for  G.  .  and  G  *  In  the  corrputer  program. 
Substitution  of  Eq.  (D-8)  into  (D-b)  yields 


A  .  1 

Aa  -  ^ 


1G21* 


"°12  G22* 


■  A6  =  D 


V  -°11 

G32*  -°12 


where 


D  = 


G  *  G  * 
u31   21 

G  *  G  « 
U32  22 


The  break  frequency  and  damping  factor  were  corrected  by  adding  one  h 
of  the  calculated  correction  factors.  These  values  of  a  and  6  were  then  used 
to  calculate  new  values  of  A  and  B.  The  process  was  continued  until  the  itera- 
tion criterion  was  satisfied. 

These  equations  were  programmed  in  FORTRAN  IV  language.  The  sub-programs 
which  perform  the  Newton-Raphson  iteration  were  originally  written  by  Mingle 
(13).  Minor  modifications  were  made  by     author  in  order  to  apply  the  sub- 
programs to  this  particular  problem.  The  source  program  is  listed  and  the 
logic      im  is  shown  in  this  appendix.  Table  D-l  defines  the  symbols  used 
in  the  program. 

The  input  to  this  program  consisted  of  the  output  of  the  program  descri- 
bed in  Appendix  C.  One  additional  card,  the  initial  guesses  for  a  and  6,  was 
read  in  as  the  last  card  in  the  FORMAT  (2F10.*0.  The  iteration  accuracy  cri- 
terion was  included  in  the  program  itself  and  was  not  ira- 
meter.  The  output  of  the  program  list:        itlon  results  thr<        ;n- 
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line  printer.  The  first  line  of  output  lists  NUM,  NPTST,  RDT,  GAMMA,  VAR, 
and  AVE  which  identify  the  output  with  the  experimental  data.  The  second 
line  of  output  gives  the  first  estimates  of  the  parameters  A  and  B,  followed 
by  the  first  estimates  of  a  and  8  in  the  third  line.  The  estimates  of  the 
parameters  are  printed  out  after  each  iteration  until  the  iteration  criterion 
has  been  satisfied.  The  last  few  lines  are  the  frequencies,  the  experimental 
power  spectrum,  and  the  corresponding  calculated  power  spectrum  in  tabular 
form.  The  final  line  states  the  squared  error  associated  with  the  least 
squares  calculation. 

The  program  took  approximately  forty  seconds  per  iteration  for  a  set  of 
sixteen  data  points.  The  number  of  iterations  depended  upon  the  accuracy  of 
the  initial  guess.  In  general  it  took  on  the  order  of  forty  iterations  to 
converge.  The  slow  convergence  was  due  to  the  changes  made  in  A  and  B  during 
each  iteration.  The  program  tended  to  diverge  if  the  initial  guess  was  too 
far  from  the  correct  value.  Appropriate  upper  and  lower  limits  were  chosen 
for  a  and  6  and  the  program  was  automatically  terminated  if  these  values  were 
exceeded. 


Table  D-I:  Input  data  and  variables  required  for 
the  IBM-1^10  break  frequency  and 
damping  factor  computer  program 


Symbol  Explanation 


A       Correlated  noise  constant  of  the  power  spectrum 

AL      First  guess  of  the  break  frequency 

B       White  noise  constant  of  the  power  spectrum 


Table  D-I  (contlr. u 


Symbol  Ex[  Sanation 


BE  First  guess  of  the  dam;     ictor 

5  Number  of  data  cards  to  be  re 

0  Frequencies  at  which  the  power  spectrum  has  been  computed 

S  Input  values  of  the  power  spectrum 

SC  Calculated  power  spectrum  using  least  squares  value  of  A,  b,  a,  and  6 

T  Statement  function 

XAL  Least  squares  value  of  the  break  frequency 

XBE  Least  squares  value  of  the  damping  factor 

NUM,  NPTS,  ROT,  GAMMA,  VAR,  AVE  (These  variables  are  not  used  directly  in 

this  program.  They  are  computed  in  the  power  spectrum  program  and  are 
used  here  as  run  indices  only.) 

SUBROUTINE  FIT 

C      Matrix  elements  used  in  calculation  of  A  and  B 

DM     Iteration  correction  factor  for  A 

GEE    Iteration  correction  factor  for  B 

Iteration  correction  factor  for  AL 
DZ     Iteration  correction  factor  for  BE 
G      Matrix  elements  used  in  calculating  DB  and  DZ 
W      Weighting  factors  used  in  the  calculation  of  A  and  B 


Pig.  26:  Logic  diagram  for  break  frequency 

and  damping  factor  computer  program. 


Start : 
Read  in  program 
and  constants 


Head : 
Data 


Output : 
0,S,SC,SUM 


Call: 
FIT 


Calculate 

SC 

SUM 


Logic  diagram  of  subroutine  FIT 


Accept 
Argument: 


Write : 
A,B 


Calculate: 
DB,DZ 


Calculate 


Calculate 
DAA,Dhb 


Calculate: 
C 


Calculate 
A,B 


1 

Set : 
XX(1)*BE 

XX(2)=AL 

Call: 
FUMCT(XX(1)) 
FUNCT(1.0bXX(l)) 

Set: 
XX(1)=BF 

'l 

Calculate: 
G 

Set : 
XX(2)=AL 

Call : 

FUNCT(1.05XX(2)) 

Logic  diagram  for  SUBROUTINE  FUNCT 


Accept 
Arguments 


EJ 


• 



.... 

I 

.  •  i   »1 \  .2.5X,'  .  ,A»5X»F10.7) 

.   i  K»I4»8X.F4.0i   •   .  •  ( i    ...  ( t F 1  . ^ ) 
•,...,  ]  3HFXPER]  •       . 
.  1 7  X  i  "  .  .  <  i  |  l  .  9  ) 
(i 
' .  i         .       »W#( 2»*B*B-1« J ►W** A 

9     '    '  .  i '   'c  »RD1  •   '    •  .VAR  .AVT 

.  03)NUM»NPTS»RDT»GAMMA»VAP,AVC 

^0(1»101)  (  (C(  I  )  »S<  I  )  )  ,  1  =  1  ,NC! 
>  A  L  •  B  E 

■  '..  L  F I T  (  X  AL  t  > 
SUM-  . 

:  i  i«itNc 
r  c  i  i  \-  ►  >   » c  <  j  n 

]  *  l+(S(I)-$CU))«*2 

:     .  (*) 

.205) (  (C( I  ) »S(  I  )  »SC(  I  )  )  *  1  =  1  >NC( 
i   6  )  SUVl 

STOP 

\0 


F  F I T  ( X*  » >  r  ) 

c    thic       r 1 1  1      TCN-RAPHSCN  correction  FACTORS 

r~r  -     rA<  fpfoutncy  and  da'-'Pint,  factor 
[25)tC(25)»G(3»2)»XX(2l  .  C ( 5  )  » W ( ?  c ) 
StOtAt  •  -»D,E»NCDS 
•  T(  10X.2E14.8) 
20?  F     '(  15X.2E14.8) 

T(   i  •  •    "  4  +  2.*A*A*Z*Z*( 2.*B*B~1 .  )+Z**4 

I C  1  =  2.' 

yd  )«i./S(  I) 

"I 

OLDB=« 

1  1=1.5 
1  C(  I  )=  . 

:  2  I-  iNCOS 
TERM- 1 • /T ( AL i  B  »C ( I  > I 
[  1 )  •  C  ( 1 )  +  T  £   •  '  I  ) 

.  (  I  ) 
»W(  I  ] 


100 


41  1 


300 


3n-i 


C  (5 
C(4 

A=( 

B=( 

DAA 

DRR 

WP  I 

XX  ( 

XX  ( 

DC 

CAL 

XX( 

XX( 

XX( 

CAL 

XX( 

DO 

DO 

G(J 

DZ  = 

DB  = 

AL  = 

IF( 

6F  = 

IF( 

WRI 

IF( 

IF( 

IF( 

IF( 

XAL 

XF>P 

P^T 

FMD 


)=C 
)=C 
C(l 
C<3 
=  A- 
=  p_ 

TP( 
1  )  = 
2)  = 
300 
L  F 
1  )  = 

1  )  = 
2)  = 
L  F 

2  )  = 
3^1 
3^1 
♦  I  ) 
(-G 
(-G 
AL  + 
ARS 
BE  + 
ABS 
TF( 
ARS 
ABS 
ARS 
A  B  S 
=  AL 
=BE 
URN 


( 5 ) +W( 
(A) +S( 
)*C(4) 
) *C ( 1 ) 

:lda 

DLDB 

3*201) 

RF 

AL 

1  =  1,2 
UNCTIX 
XX(  1  )* 
BE 

XX( 2  )* 
UNCTIX 
AL 
J=2,3 
1  =  1,2 
=(G< J, 
(1,1)* 
(1,1)- 
DB*0.5 
(AL)  .G 
DZ*0.5 
(BE  )  .G 
3,202) 
(DB/AL 
(DZ/EE 
(PA  A /A 
(PBB/B 


I  ) 

I ) *W ( I ) 

-C(3)*C|5))/IC{1)*C(1)-C|5)*C(2)) 
-C(4)*C(2  )  )/(C{ 1 )*C( 1 )-C(5)*C(2 )  ) 


A,R 


X(2 ) ,XX( 1 ) ,G(  I  ,1  )  ,G( I  ,2)  ) 

1.05 

1.05 

X(2) ,XX( 1 ) ,G(3,1 )  ,G(3,2)  ) 


I )-G( 1  , I  )  ) /( .05*XX( J-l  )  ) 

G(3,2  )+G( 1  ,2)*G(3»1  )  )/ (G(?,l  )*G(3»2)-G( 2,2)*G(3»1  )  ) 

G(2,l )*D7)/G(3,]  ) 

T.5. )RETURN 

T.2. JRETURN 

AL,BE 

) .GT.l .E-3)GCTC499 

) .GT.l .F-3)GCTC4P9 
) .GT.l .E-3)GCTC499 

) .GT.l .F-3  )GCT0499 


SUB 
THIS 
D  I  M 
COM 
T(A 
U(  A 
V(  A 
Pl  = 

P2  = 
PC 
Rl  = 
R2  = 
RFT 
FMD 


ROUT 

PRO 

ENS  I 

MCN 
,  R  ,  W 
,  B  ,  W 
♦  R  ,  W 
0. 
C. 

1  1  = 
Rl 
R2 
URN 


INE  FUNCT( AL,BE,R1,R2) 

GRAM  COMPUTES  THE  VALUES  OF  THE  NORMAL  EQUATIONS 

ON  S( 25), 0(25) 

S,C,A,R,X,Y ,NCDS 

)=A**4+2.*A*A*W*W*(2.*B*B-1.)+W**4 

|s4«*A**3+4«*A*W*W*(2«*B*B-l«  J 

)«8«*B*A*A*W*W 


2,NCDS 

+(S(I)-B-A/T(AL,BE,C(I)))*U(AL,BE,C(I))/(T(AL,RF»C(I)))**? 
+(S( I )-B-A/T( AL,BE,C( I ) ) )*V (AL,BE,C( I ) ) /( T( AL,BE ,C( I ) ) ) 
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APPEf. 


•  •  :     N    FUNCT ' 
fRUM       VLCULATICf 


'•-  " 
F II  "24. 

.    DATA    POINTS    =    1986 

Ml NG    CCEFF     =    0.0 

[GHT  ■  6.346 
VARIANCE  ■  »1497 


Tivfr    LAG 


rc  fx\< 


FREQUENCY 


spfctpu" 


.Oo 

.      -166 

•  oe 

.16666 
.20 

. 
.29166 

.416^6 


oooor 

.  '0106219 
.72083944 
.  1274946 
•31102967 

.  14^8*018 
.02417163 

-.05000729 

-.0  750660/. 

-.  -  )99005 

- .    283<J : 
•02127231 
•  ■-.  c  5 


.0000 

1  .r 

2*0000 

^.0000 
4.0000 
.  ^000 
6.0000 
7.0000 
8.0000 
9.0000 
10.0000 

1 1 .oooo 

12.oono 


.211948R40 
.20081 1200 
.  1  18424960 
.0400454 

.010031101 
.005036702 
.002628733 

.001969492 
.001572291 
. 00124372^ 
.000924842 
,OO0PQ-*127 
.000887  7  64 


APPENDIX  E»  CONTINUED 


RUN  NUMBER  A-8 

FILM  SPEED  =  32. 

NO.  DATA  POINTS  =  2020 

PREWH ITEMING  COEFF  =  0.0 

AVERAGE  RED  HEIGHT  =  6.3! 
VARIANCE  =  .1628 
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TIME  LAG 

AUTO  FXN 

.00000 

1.00000000 

.03125 

.84298746 

.06250 

.55256338 

.09375 

.25027385 

. 12500 

.01394600 

.15625 

-.12175533 

.18750 

-.  14894121 

.21875 

-.09048296 

.25000 

.01772836 

.28125 

.  11886211 

.31250 

.16^63701 

.34375 

.  15925221 

.37500 

.09887175 

.40625 

.01577766 

.43750 

-.05423553 

.46875 

-.09633460 

.50000 

-.10895160 

FREQUENCY 

SPECTRUM 

.0000 

.125815070 

1.0000 

.115548200 

2.0000 

.111778880 

3.0000 

.105264380 

4.0000 

.059169837 

5.0000 

.017805296 

6.0000 

.008420228 

7.0000 

.005805843 

8.0000 

.003524270 

9.0000 

.002274329 

10.0000 

.001699633 

11.0000 

.001459364 

12.0000 

.001140298 

13.0000 

.000926287 

14.0000 

.000885206 

15.0000 

.000916296 

16.0000 

.000948224 

APP  ND1X  £ i 


•  -10 

■  ■•  . 


■  •  FREQUEf  SPECTRUM 

.  000  .  I6?0  7r-  • 

.02083       .9361  1 .oooo  . 1 

.  -.166        .82872840  2.0000  .128335370 

.  •  almoin  - .   no  .078231  c'07 

.  3148111  >.0000  .028597595 

.  7Q761  P  1  ^.0000 

.  1  .        .'t1  6.0000  »0(      253 

.14!  . '  7.onoo  .003! " 

P. 0000  .00214767^ 

.•75"        -.'•  .  .001375127 

10.0000  .001082971 

16        -.  ll.'H.OO  .00080F 

12.0000  .000654?!^ 

.   "  B3       -.06854996  .0000  .000551840 

L66       -.04052897  14.0000  .000563316 

.  1250       -.00786  15.0000  .000574082 

.  •  .  .00049^ 

.^5416         .        ^7  17.0000  .000426" 

. 375  .06695139  IP .0000  .    430470 

lo.npoo  .000434564 

.  .  9149A  20.^ooo  ,oooA?QtU 

,09180393  21.0<  .000480051 

.  B3612'  .   00  .000557] 

.47916        .0<    >67  23.0000  .000525233 

.5000C         .   >24193  24.  .000469041 


APPENDIX  E,  CONTINUED 


RUM  NUMBER  A-ll 

FILM  SPEED  =  48. 

NO.  DATA  POINTS  =  2  00  7 

PREWHITENING  COEFE  =  0,0 

AVERAGE  BED  HEIGHT  =  5.473 

VARIANCE  =  .0624 
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TIME  LAG 


AUTO  FXN 


.ooooo 

i.oooooooo 

.02083 

.88257535 

•04166 

.74453082 

.06250 

.58269833 

.08333 

.41632020 

.10416 

.26263229 

.12500 

.  13111323 

.14583 

.02413762 

.16666 

-.05375740 

.18750 

-.09755576 

.20833 

-.11250285 

.22916 

-.0Q327403 

.25000 

-.06366697 

.27083 

-.03092280 

.29166 

•00238947 

.31250 

.02931^61 

.33333 

•058  264  7  2 

.35416 

.08335522 

.37500 

.09966050 

.39583 

.09444705 

.41666 

.08535229 

.43750 

.06823623 

.45833 

.05392727 

.47916 

.0367676^ 

.50000 

.01725314 

REQUENCY 

SPECTRUM 

.0000 

.134356380 

1.0000 

.128890900 

2.0000 

.118^-77220 

3.0000 

.086477674 

4.0000 

.039560181 

5.0000 

.014869063 

6.0000 

•008029430 

7.0000 

.005981229 

8.0000 

.005044259 

9.0000 

•003222917 

10,0000 

.002261517 

11  .0000 

.002165^60 

12.0000 

•0H196737Q 

13.0000 

•0O]810653 

14.0000 

.00162468* 

15.0000 

.001456877 

16.0000 

.00154^6^8 

17.0000 

•001519272 

18.0000 

.001298242 

19.0000 

.001149999 

20.0000 

.001179872 

21.0000 

.001323238 

22.0000 

.001267264 

2^.0000 

.001132204 

24.0000 

.001136224 

APPf-MMX   . 
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.  29743  11.0"  .  00666<  • 

.0034  12.0(  ,0r     H89 

702268  13.0000  .000580676 

06       -.00326367  14.0000  .  ^0565 

.  1250        .  1150306  15.^000  .000509916 

.   ••.-6175  16.0000  .  ^0421011 

.  54  16        .06727^  17.0000  .000391- 

=c  18.nr  .000420968 

lQ.oooo  ^nnoupo^-'P 

20.0000  .000543^1^ 

.  3272'.  21.0000  .000521 

33        .198412  2?.0oon  .    00424! 

.0000  .  '  ^50432 

.  ">0984-  24.0000  .000338753 
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Addpmdix 


CONTINUED 


RUN  NUMBER  R-l 

FIL'<4  SPEED  =  48. 

NO.  DATA  POINTS  =  199] 

PREWHITENING  CC^FF  =  0.0 

AVERAGE  BED  HEIGHT  =  6.782 

VARIANCE  =  .0265 


TIME  LAG 

AUTO  FXN 

.00000 

1,00000000 

.02083 

.POP 916 8 6 

.04166 

.77911783 

.06250 

.63394203 

• 08333 

.48198448 

.  1 C  A  ]  6 

. 34584076 

. 1250O 

.22656614 

.14583 

.12667270 

.16666 

.05777568 

.18750 

. 01582368 

.20833 

-.001113  15 

.22916 

.00896590 

.  2  5  0  C  0 

.040448  8  7 

•27083 

.  "80104  5  3 

.291 66 

.13014352 

.31250 

.  1727Q966 

.33333 

.20^75190 

.  3  5  4  1  6 

.22457987 

.37500 

.23444710 

.39583 

.23226245 

.41666 

.22025802 

.43750 

.20495119 

.45833 

.17278218 

.47916 

.14087727 

.50000 

.  11195771 

FREQUENCY 

1.0000 
2.0000 
3.0000 

4.0000 

c.0000 

6.0000 

7.0000 

8.0000 

9.0000 

10.0CC0 

11.0000 

12.0000 

13.0000 

14.0000 

15.0000 

16.0000 

17.0000 

18.0000 

19.0000 

20.0000 

21 .0000 

22.0000 

23.0000 

24.0000 


SPECTRUM 

. 17587031  o 
•1414376 
. 106494440 
•079748284 

.P"^9P4464 
.011988531 
.007950770 
.005410273 
.003585455 
.002646984 
.002169852 
.001699634 
.001474675 
.001291941 
.001162533 
.001230864 
.001430585 
.001 386754 
.001223664 
.001093622 
.001101756 
.001056239 
.000942618 
.001007170 
.0010880^6 


•    •  IX    F  . 


■    1983 

,0 
7  .o  19 
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\G 

MJTC 

. 

1  . 

. 

.  P402U 

•  o/»: 

.82253: 

:5C 

.71233182 

. 

.   9647]  Rfl 

.10.. 

•  • 

• 

. 4044??Q6 

. 

. 

• 

•  2761 

.09 

. 

.  '  ' 

. 

. 20916504 

•  250 

.  ' 

83 

.23186- 

. 

»65 

.  ' 

3  7  1 

. 

JO  4  9 

.35 

.: 

. 

•  7? 

. 

. 

• 

.  ' 

. 

.  0449785 

• 

. 1 7325690 

16 

.  1426K 

. 

.11002082 

FREQUEI 

00 

1 . 

2.0000 

3.0000 

^.0000 

.^000 

.0000 
Q.0000 

1 0  •  oc 

11 .0000 

.00 

13.0000 

1 4.0000 

.0000 
16.0000 
17.0000 
18.00 
19. 

.   '00 

00 

.0000 

24.0000 


.   1   ■• 

.2^     "90 
. 169350500 

.086123 
.056621629 

•  »62 

.010260: 

■  1  ' 

•  r86 

.  "22^270? 

.        15671- 
.     "13  7  2113 
.001  116717 
.00101 2 16^ 
.     : 10786  " 

.noi  01  p.. 

.  n  n  n  o  6  0  <4  3  ? 

.0009*7 
.nnoo/47> 

,0008417 
•00067 

.  !  0  6  2 c  : 

'08 


APPENDIX  E,  CONTINUED 


PMM  MuMfipp  q_^ 

FILM  SPFFD  =  ^fi. 
NO.  DATA  POINT?  =  1984 
PREWH ITEMING  CCEFF  =  0.0 
AVERAGE  BED  HEIGHT  =  9.096 
VARIANCE  =  .0581 
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T^E  LAG 

AUTO  FX\ 

.CO  000 

1.00000000 

.  0  ?  0  8  3 

.89834163 

.04  166 

.80451764 

.06250 

.69560239 

.08333 

.57725250 

.10416 

.46384292 

.12500 

.35471477 

.14583 

.25752721 

.16666 

.18204503 

.18750 

.12730112 

.20833 

.08886573 

.22916 

.C6694679 

.25000 

.06697204 

.27083 

.07969779 

.29  166 

.  10366097 

•31250 

.118  8  814  2 

.13333 

.13710367 

.^5416 

.14655181 

.37500 

. 14982124 

.39583 

. 14069492 

.41666 

.12178577 

.43750 

.10221831 

.45833 

.03300623 

.47916 

.05797830 

.5000,') 

.0309567*5 

FREQUENCY 

SPECTRUM 

.^000 

.203897860 

1.0000 

.163518340 

2.0000 

.  103799630 

3.0000 

.060774819 

4.0000 

.023080472 

5.0000 

.008225664 

6.0000 

.006110788 

7.0000 

.004220359 

8.0000 

.003131932 

9.0000 

.002665458 

10.0000 

.002633468 

11 .0000 

.002250778 

12.0000 

.00] 656465 

13.0000 

•0O1489545 

14 .0000 

.001 541648 

1 5.0000 

.001544183 

16.0000 

.001 5  45  741 

17.0000 

.001538828 

18.0000 

.001437685 

19.0000 

.001309597 

20.0000 

.001360221 

21.0000 

.001305909 

22.0000 

.001105277 

23.0000 

.001  156490 

24.0000 

.001295436 

Apr 


:  • 


',6 


A 1 1  T  3     ■ 


;  NCY 


.        J83 

.16666 
.22 

.      :83 

416 

.        583 
•  4 1  <  • 


.     6084C 
.7  .     77 

1  05 
'1251 

-.0063K 
-.  "£  **7949 

-  .      1 4  ?  7  6  8 1 

-  .     7  7  7 1 7  7  7 
.  >3827 

-•04 

-.  016491 
•  c 

•  0  1 

•  2° 
56 

187 
•  < 

.    ••     >5012 

.     2678! 


.     0  n  0 

1.0000 

.0000 

.     000 

4.0000 

5.0000 

9.0000 
10.0000 
11.0000 
12.0000 
1  '-.0000 
14.0000 
15. 

16.0000 
17*0000 
1 8  •  Of 
19.0000 

.   ^00 
2?. 

23.0000 
24.0000 


. 
•04624^398 

.002371478 
.  -21181 
.   2014084 
.   Ib61 
. 

71 
. 

•00173 
.     'U4434 
.     114084 

•  0014 
.001465' 
.     51537] 
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APPENDIX  E»  CONTINUED 


PUN  NUMBFR  R-5 

FILV  SPFFD  =  32. 

MO.  DATA  POINTS  =  2027 

PRFV.'H  ITEMING  COFFF  =  0.0 

AVERAGE  BED  HEIGHT  =  6.272 

VARIANCE  =  .0436 


TIME  LAG  AUTO  FXN  FREQUENCY  SPECTRUM 

.00000  1.00000000  .0000  .214457530 

.03125  .87946292  1.0000  .172492240 

.06250  .71687230  2.0000  .105259390 

.09375  .5469°361  3.0000  .055077043 

.12500  .39343991  4.0000  .02158006^ 

.15625  .26407762  5.0000  .010305236 

.18750  .16543217  6.0000  .006361125 

.21875  .10691922  7.0000  .004666722 

.25000  .07837673  8.0000  .003696520 

.28125  .07196643  9.0000  .002644039 

.31250  .07754425  10.0000  .002010750 

.34375  .09801685  11.0000  .001806813 

.37500  .12095030  12.0000  .001626639 

.4C625  .1371652'3  13.0000  .001576614 

.43750  .140^586  14.0000  .0^1571644 

.46875  .15791032  15.0000  .0014312^1 

.50000  .15166118  16.0000  .001330199 


APP  • 


.0 


LAC                 AUTC    F>  FREQUENCY                      ECTRl 

•  Of 

•o;                     .  .00 

.04166  .8101 

•  4838791  U  A.or-oo  .o  • 

.    '     ^m?3  5.0(  op 

.  \6Ut  .0^00  .                  7^'. 

•           33                      .  7.0000  .     •   - 

30  .00218? 

-.117924  .0000  .001469301 

-.1502437^  10.0000  .001242401 

.     2016                 -.  15200  11 ,0000 

-.l?flflll70  12.0000  »0( 

.27                         - .      • 1  439]  ]  .on  . 

-.04556423  14.  .00061 

.002671  lc  .^r^nr\  .     '           1485 

•33333                   .045551  16.0000  .        )717 

,08031584  .       on  .ono^ii 

.100679^  .  .            t78865 

.i057nm  19.0000  .oooa?^ 

.      >581901  20.0000  .000430322 

.     74807!  21.0000  .     00435 

.     -4117206  22.0000  .000-      ■ 

16                     .     0146686  .00  . 

137  24. or  .0004 


APPENDIX  E>  CONTINUED 


RUN  NUMBFR  C-2 

FILM  SPEED  =  48. 

NC.  DATA  POINTS  =  1985 

PREWH ITEMING  COEFF  =  0.0 

AVERAGE  BED  HEIGHT  =  6.556 

VARIANCE  =  .1366 
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TIME  LAG 


AUTO  FXN 


,00000 

1 .00000000 

.02083 

.94812779 

.04166 

.85029046 

.06250 

.72042544 

.08333 

.57481385 

.10416 

.42788686 

.12500 

.28904537 

.14583 

.  16498^66 

.16666 

.06017877 

.18750 

-.0106018° 

.20833 

-.076468] 6 

.22916 

-.11100705 

.25000 

-.  11960208 

.27083 

-.109351 8? 

.29166 

-.08157044 

.31250 

-.04302687 

.33333 

-.00040129 

.35416 

.04130599 

.37500 

.07714757 

.39583 

.10249287 

.41666 

.  11701291 

.43750 

. 12188229 

.45833 

. 1 1502798 

.47916 

.09983006 

.50000 

•07867914 

FPFQUENCY 

SPECTRUM 

.0000 

.170286580 

1.0000 

.162035830 

2.0000 

.133311480 

3.0000 

.073915408 

4.0000 

.021932241 

5.0000 

.007773081 

6.0000 

.004792393 

7.0000 

.002433823 

8.0000 

.001572545 

9.^000 

.001131297 

10.0000 

.000896976 

11.0000 

.000659970 

12.0000 

.000498391 

13.0000 

.000425175 

14.0000 

.000434024 

15.0000 

.000448404 

16.0000 

.000420777 

17.0000 

.000306996 

18.0000 

.000264928 

19.0000 

.0H0MP856 

20.0000 

.000342534 

21.0000 

.000299158 

22.0000 

.000250441 

23.0000 

.000259987 

24.0000 

.000279^1 1 

APP         »  CONT  INI. 


a     j  <; 

FF  ■  0 • C 

■  . 
VA-  |AI      ,23! 
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IE  L 


A  U  T  C 


FRFOUFNCY 


SPECT' 


.0? 

.   -166 

.  14583 
.16666 

.    750 

•   333 

.22 

. 

.27 

166 
.3] 
.^"^ 

S66 

.47 
.500 


on 

.  '78734 

. 

•46683] 
.  76631 

•04 
71 
.185686 

5  4  ] 

.12515^ 

,134 '  • 
.1401438? 

.  6 1  n  7 1 

,  1(      ■ 9  1 

. 17740397 
. 17587724 

.  L69711 


.  ooo 

2.0000 
3.0000 
4.0000 

6.0000 

7.0000 

8.0000 

o.Or 

.   00 

1 1 .0000 

12.0000 

13.0000 

lA.r^oo 

1 5.0000 

.   on 

1  7.0000 

18.0000 

lQ.r^oo 

?o. or 

.  00 
.  .00 
.0000 


.2'-     -20 
. 185902? 
.10191111 
.  501033 
•0168] 

.Or 

^9440 
.002^6521 
•  001863;  • 

.001409158 
.000854139 
.OO069R390 
.000788422 
.   0  7  4 18 77 

mi2 

.0006221 

8? 
.        042'^ 
.     "0^07084 
.      104       U 

.  ^1401382 
.00041 1947 


APPENDIX  E,  CONTINUFD 


RUM  NUMBER  D-l 

FILM  SPEED  =  48. 

NO.  DATA  POINTS  =  2074 

PREWHITENING  CCEFF  =  0.0 

AVERAGE  BED  HEIGHT  =  8.236 

VARIANCE  =  .0491 
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TIME  LAG 

AUTO  FXM 

.00000 

1.00000000 

.02083 

.94111165 

.04166 

.83507594 

.06250 

.70498429 

.08333 

.56893305 

.10416 

.44019811 

.12500 

.32496865 

.14583 

.22903111 

.16666 

.15719094 

.18750 

.10712258 

.20833 

.07702486 

.22916 

.06939048 

.25000 

.080030  91 

.27083 

. 10074724 

.29166 

•  12480649 

.31250 

.14772051 

.33333 

. 16889458 

.35416 

.18620364 

.37500 

.19934555 

.39583 

.20425294 

.41666 

.20202880 

.43750 

. 191050S9 

.45833 

.  1  7  3  A  9  2  3  2 

.47916 

. 1 53988  8  7 

.50  000 

. 13382857 

FREQUENCY 

SPFCTRUU 

.0000 

.205496330 

1.0000 

.162367470 

2.0000 

.106652920 

3.0000 

.068475306 

4.0000 

.028053541 

5.0000 

.010340945 

6.0000 

.005919274 

7.0000 

•003844248 

8.0000 

.007442703 

9.0000 

.001462286 

10.0000 

.0^1265954 

11.0000 

.0010929R8 

12.0000 

.000807471 

13.0000 

.000550864 

14.0000 

.00045975Q 

15.0000 

.000459497 

16.0000 

.000480860 

17.0000 

.000447790 

18.0000 

.000392917 

19.0000 

.0003419Q8 

20.0000 

.000310482 

21 .0000 

.000324661 

22.0000 

.00030^31  3 

23.0000 

.0002919?? 

24.0000 

.000303261 

APPFNDIX  E»  CHN7 INUf 0 


•   '   '  <FR  F-1 
F II        ">  ■  4  8  . 

PREWH] TUNING  fCFFF  =  0. 

T  ■  10.  • 
•  E  -  .051 7 


TIME  LAG 

AUTC  FXN 

. 

1 .00000000 

. 

.95100774 

.  4166 

.87043467 

»0< 

. 7683^741 

. 

.  ^>  5  2  3  6 1 

.  '   16 

.  ^38241  ^8 

. :2500 

.43523666 

.14583 

.34296695 

.16666 

.26954763 

.  18750 

.22151888 

.208 

.  19390578 

.22916 

.  18866029 

.2! 

.204504  2  2 

B  3 

.23095003 

.29166 

.26507^1 l 

.31 

.  ^"0414068 

,33333 

.  13739869 

•  -*M6 

.  \6  1  -5^n«^ 

•  375 

•37110104 

.3931 

.^6S74207 

.41666 

.  )4 992891 

.43750 

.32651981 

.45833 

.29503550 

.  -7916 

. .6684460 

30 

.24235140 

FREQUENCY 

.0000 
1  .or.no 
2.0000 
■  .onoo 
4.0000 
5  .r^ooo 

^.0000 
7.0000 
8.0000 
9.0000 
10.0000 
11 .0000 
12.0000 
13.0000 
14.0OQ0 
15.0000 
16.0000 
17.0000 
18.0000 
19.0000 
20.0000 
21.0000 
22.0000 
23.0000 
24.0000 


SPECTRUM 

•252015390 
.  176176300 

.093446826 
.06156**92 
•0211321 
•005585339 

.OO^02172P 
.002^084^0 

.00]  702"' 

.  91165651 

.000791363 
.000652180 
.000587 7 
.000489721 
.OA054O? 

.0O065396R 
.00065842 1 
•000458661 

.000324691 
.000329287 
.000330484 
.000313- 
.00031 1  • 
.  ^0301892 


APPFNDIX  Et  CONTINUFD 


RUN  NUMRFR  A-10-R 

FILM  SPEED  =  48. 

NO.  DATA  POINTS  =  1987 

PREWHITENING  CCEFF  =  0.0 

AVERAGE  BED  HEIGHT  =  6.461 

VARIANCE  =  .1453 
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TIMF  LAG 

AUTO  FXN 

.00000 

1 .00000^00 

.02083 

.92594445 

.04166 

.82217585 

.06250 

.70320785 

.08333 

.58109976 

.10416 

.46072551 

.12500 

.34639267 

.14583 

.24368361 

.16666 

.  15441976 

.18750 

.08804333 

.20833 

.04152210 

.22916 

.02004466 

.25000 

.01571293 

.27083 

.02220691 

.29166 

.03095549 

.31250 

.04371305 

.33333 

.05503687 

.35416 

.06518827 

.37500 

.07463563 

.39583 

.07994730 

.41666 

.08200350 

.43750 

.07^82060 

.45833 

.07707245 

.47016 

.07414934 

.50000 

.07259227 

FREQUENCY 

SPECTRUM 

.0000 

.194418960 

1.0000 

.166139080 

2.0000 

.113465500 

3.0000 

.062388525 

4.0000 

.023223737 

5.0000 

.008370066 

6.0000 

.005547935 

7.0000 

.004704718 

8.0000 

.003548422 

9.0000 

.002375362 

10.0000 

.001812077 

11 .0000 

.001545485 

12.0000 

.001308950 

13.0000 

.001063812 

14.0000 

.000915109 

15.0000 

.000884415 

16.0000 

.000808838 

17.0000 

.000678137 

18.0000 

.000622977 

19.0000 

.000621921 

20.0000 

.000625028 

21 .0000 

.000675774 

22.0000 

.00067227^ 

2^.0000 

.000552688 

24.0000 

.000479305 

IUI  l' 


1  -p 

0.0 
RED    HFIGH1     -     • . 
V*  ■    .031 


AUTC 


FkFUUf-NCY 


•  00' 

|  66 

.  •      416 
.  1  ?  ' 

.      =.e? 

. 16666 

•  22 
.2b. 

.?9166 

.-Al 

.  3^416 

583 
.4  1666 

.45833 

.47 

100 


1  .  >00 

.  •    >240 

. 78? 3  "4  80 
.  1  ?1 06  4? 
.  •  ^718 
.     ■ 

.  1  1219224 
77200 

-. 1 5279307 
-.  L69 

-.  lr.?28250 
-.  1260?' 
-.07258031 
-.01597424 

.     • .  •  0  8  7  0  U 

•  131  M687 
.  1  5743525 

.  16316701 
.  15377730 
,130 

•  0    5657 

.  5209735 

. 


.  00 
1 .0000 
2.0000 
3.rono 
4.0000 
«^.oooo 

6.0000 

7.0000 

8.0000 

9.0000 

10.0000 

1 1 .0000 

12.0000 

13.0000 

I4.oono 

l r .0000 

16.0  000 

17.0000 

.Of, 00 

19.0000 

20.0000 

.0000 

22.0000 

23.0000 

.   00 


. 1 30466020 
. 12902045C 
•129344 

•03821  • 

.  0  1  ?  6  c  7  3  7  Q 

.003014030 

. 

.0< 

,001  1  7<  ^65 

•  OOOS 
.0008151  4 1 

•000903852 

.or 

,ooo>-. 

•  00061 

.00077 1   - 

•  0  0  0  fl 

.000654 <81 

•  000^ 

•Of     '  '71 

] 
•0007- 


APPENDIX  E»  CONTINUED 


RUN  NUMBER  C-2-R 

FILM  SPEED  =  48. 

NO.  DATA  POINTS  =  1990 

PREWHITENING  CCEEF  =  0.0 

AVERAGE  BED  HEIGHT  =  6.618 

VARIANCE  =  .1670 
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TIME  LAG 

.00000 
.02083 
.04166 
.06250 
.08333 
•10416 
. 12500 
.14583 
. 16666 
. 1 8750 
.20833 
.22916 
.25  000 
.27083 
.29166 
.31250 
.33333 
.35416 
.37500 
.39583 
.41666 
.43750 
.458^3 
.4 70  16 
.50000 


AUTC  FXN 

1.00000000 

.93264920 

.82983190 

.69472834 

.54832641 

.40417730 

.26704779 

.14778225 

.052^0020 

-.01535609 

-.05563374 

-.07068705 

-.06  54  5547 

-.04?99218 

-.01437944 

.01681533 

.04608173 

.06983401 

.08910532 

.1027270^ 

. 1 1074633 

. 10940901 

. 1000^676 

.08M0980 

.05938738 


FREQUENCY 

SPECTRUM 

.0000 

.170751920 

1 .0000 

.156231000 

2.0000 

.125227890 

3.0000 

.075791564 

4.0000 

.026885750 

5.0000 

.008333803 

6.0000 

.004758808 

7.0000 

.003329700 

8.0000 

.002306406 

9.0000 

.001686078 

10.0000 

.001354O0S 

11 .0000 

.000997002 

12.0000 

.000810012 

13.0000 

.000684352 

14.0000 

.000629196 

15.0000 

.000582746 

16.0000 

.000579399 

17.0000 

.000620113 

18.0000 

.000650301 

19.0000 

.000624279 

20.0000 

.00058968? 

21  .0000 

.000580496 

22.0000 

.00056082? 

21.0000 

.000537960 

24.0000 

.000541696 

APl'FNDI  ■ 


■ o g  c  KIT 


•   •        ■  -7 

■  1  • 
DAMPING  FACTOR  ■  ,6232 
ERROR  SQUARED  =  .00002169 


ENCY         EXPERIMENTAL  CALCULATFD 

1  .                .?  10H1  1290  ,?on^7^?o 

•60  . 1)91 68^30 

.  ooo            ,o/<o^/,^647  .o     ^389 

4.  ^00            .o1    •  ]  oi  ,0l240f 

,0(    •  702  .0054 

.   00             .002628733  .00295^76? 
7»0000            .00  1969492 

8.0000            .001572291  .001382346 

9. mOO            .001243723  .001112292 

.000924842  .P009S8790 

11.                 .on  nfl  03 127  .       •  n4A 

I   . '00'               .000887264  ,000807984 


■   ■■'PR  A-R 
BREAK  FREQUENCY  =  1.235 

DING  FACTOR  =  .c184 
ERROR  SQUARED  =  .000314"-7 


FREQUENCY  EXPERIMENTAL  CALCULATED 

l.OOC  .11SS4820P  .105^4^0^0 

.000  .111778880  .122285850 

►380  • 1027c6460 

►837  .050746909 

00  .017PO5296  .022037774 

00  • 00 H 4? 02 28  .0  10*78  8  04 

jO  .  '5805843  .005771 

.0000  .0035  .003527850 

.0000  .      ••  ?29  .002 

:  .  .  .  0 ' 

00  .  .00135894 

'  .  .001 1 40298 

,0000^6287  .000-   ■ 

•06  .0008 

'  •  .  0091 62^6  .        107 

.   0048224 
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APPENDIX  F»  CONTINUFD 


RUN  NUMBER  A-10 
BREAK  FREQUENCY  =  2.666 
DAMPING  FACTOR  =  .6284 
ERROR  SQUARED  =  .00011606 


FREQUENCY  EXPERIMENTAL  CALCULATED 

1.0000  .153473460  .150636840 

2.0000  .1283^5370  .134023860 

3.0000  ,078231507  .070051189 

4.0000  .028597595  .028514^64 

5.0000  .009606935  .012455075 

6.0000  .005428253  .006200226 

7.0000  .003512119  .003476474 

8.0000  .002147673  .002155951 

9.0000  .001375127  .001456802 

10.0000  .001082971  .001059527 

U-0000  .000808302  .000820501 

12.0000  .000654213  .000669779 

13.0000  .000551840  .000570951 

14.0000  .000563316  .000503974 

15.0000  .000574082  .000457283 

16.0000  .000499331  .000423929 


RUN  NUMBER  A-ll 
BREAK  FRFQUENCY  =  3.005 
DAMPING  FACTOR  =  .6084 
F^ROR  SQUARFD  =  .00007918 


FREQUENCY  FXPERIMENTAL  CALCULATED 

1.0000  .128890900  .125436160 

2.0000  .118477220  .123984830 

3.0000  .086477674  .081670601 

4.0000  .039560181  .038199782 

5.0000  .014869063  .01779366* 

6.0000  .008029430  .00^39351* 

7.0000  .0^5981229  .00*689805 

8. noon  .00*044259  .00^9^177 

.0000  .003222917  ,00294436] 

10.0000  ,00??MM7  .002406869 

11.0000  .002165960  .00208440* 

12.00  .001«67379  .00188 157* 
•   00  .001810653  .001  /  i 

]/t-  .00  K  ■  •  .0016  590  7? 

15.001  .Oi  B77  ,001  l     • 
16.  unco  •001543< 


APPENf 


•'•  -  1  2 
•  •  • 
IPING  FA  • 


•  00    I 


FRF' 


IENTAL 


CA| 


7.  J 

.  oco 

.  coo 

12.0 

.   00 

.0000 
.0000 


.1H4619630 
•120057850 

. 017909701 
.004358590 
.002322656 
.on] 707799 
•001409R85 
.OOK06859 
.^0780414 
•0006*6610 
.00061 oe89 
.OO05R0676 
.000565903 
.000609916 
.000421011 


. 183565230 
•123524240 

.043  176, 

.  0  1  5  1  <v  *  3  7  c 

.006493905 

.003302666 

.001935QPO 

.001 276C70 

.000976'. 

.000727^ 

.000607- 

.0005  "32186 

.O0O4P2508 

,0004486  . 

.000425282 

.000408466 


o :^P  o -  1 

ppr • .  :  RFQUFNCY 
FACTOR 

Fpr       )MAPFP  ■ 


i  7  .  «>  ?  fl 
.6C 
00076'179 


FRhOUf 


FXPFR  If-'ENTAL 


CALCULATED 


1 .  OOCO 
.  OCO 

1 1 .uOOO 
.  *000 


.  14  14?  7(^90 
. 106494440 

.079748284 
,0??oh^64 
•Oil  153] 
.^79*0770 
►73 

35854 
,on?646 
.007169852 

L  6  9  9  6  3  4 
. 

1  162! 

I  64 

66 


.  136792490 

•  1165 

.068629' 
.   '24Q60?4 

•015727685 

•  00844 
.On'  118^11 

•002534834 

.  oo?r . 

.001696238 
141 
.001 362; 
.001 2724 
.0012094 
.001 164^47 
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APPfNDIX  F  *  CONTINUED 


ri.im  K)\  |MR|  '?  P-2 
BREAK  FRFQUFNOY  =  2.8  24 
DAMPING  FACTOR  =  1  .01 97 
ERROR  SQUARED  =  .00016310 


FREQUENCY  EXPERIMENTAL  CALCULATED 

1.0000  .169350500  .167973660 

2.0000  .086123340  •0928824™ 

3.0000  .056621629  .046386576 

4.0000  ,025448462  .023693398 

5.0000  .010260240  .012921703 

6.^000  .OP719104O  .00760^603 

7.0000  •004314886  .00482516? 

8.0000  .003057964  .003288630 

9.0000  .002252792  .002393657 

T-.0000  •001748281  .001847975 

11.0000  ,001567869  .001501733 

12.0000  .001372113  .001274283 

13.0000  .001116717  .001120279 

14,0000  .001012189  .001013202 

15.0000  ,001078610  .000936994 

16.0000  .0OHO8327  .000881622 


RUN  NUMBFR  R-3 
BREAK  FREQUENCY  =  2.552 
DAMPING  FACTOR  =  .7748 
ERROR  SQUARED  =  .00012011 


FREQUENCY  EXPERIMENTAL  CALCULATED 

1.0000  .163518^40  .1620^7730 

2.0000  .103799610  .108775690 

3.0000  .06077401^  .0M72119O 

4.  .02^090472  .023120714 

5.0000  .008225*64  .011473161 

6.  .006110788  .006541276 

7.0000  .   '4220359  .0 

•  00313193Z  .00  )< 

9.  i  .002665458  .0024< 

1'  .0000  .002633468  .002' 

11.0000  ,002250778  .001869848 

1  .  .001656465  .00]  ' 

13.00!  .001489!  ,00]  6  I  ] 

14.  n  ,00  |  !  /,  1648  .001"  ■ 

1  5  .  .001544183 

16.  .00154574] 


APPI  .  NTJNUFD 
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■  "JFR  P-4 
•   FPFOUFMCY 
'"PING  FACTOR  = 
ERROR  SQUARED  = 


■  3.353 

.6466 
.00002195 


FREQUENCY 


FXPFRI^FNTAL 


CALCULATED 


1  .r 

00 

4  .  000 

.  ooo 

6. 0 000 

.  000 

b .uOOO 

9.0000 

10.0000 

1 1  . 0  0 

12.0  000 

1 3.0000 

14.  -^po 

•     .  ^00 

lft.'HiOO 


•118700120 

.  l^QfUlSOO 
.08*5^7468 
.046245^98 
.022121647 
.012158953 
.007156764 
.004877323 
.003823559 
.002999049 
.002371478 
.0021 18951 
14  0  8  4 
.no t 8^ ip p7 
.00187^095 
,no?09797l 


•  1 1 6  6  8  9  C 

.  1  1-2120260 
•OF   S9513 

.02311 1043 

.012550352 

.007576237 

.005080069 

.003736697 

.002967032 

.002501927 

.002207947 

.002014929 

.001 R«4^ 

.001702730 

.001727*01 


RUN  NUMBER  B-5 
BREAK  FREQUENCY  =  2.550 
OA^ping  FACTOR  =  .8265 
ERROR  SQUARED  =  .00002940 


FPEOl'FNCY 


FXPFPI^FNTAL 


CALCULATFO 


.  000 

2.0000 

.  000 

.  '000 
6  •  0  C 

.  ooo 

.   GO 
16.0000 


.172492240 
. .  5259390 

.055077043 
1580 

.  10305236 
.006361125 
.0046667/2 

1626639 

.    13  30199 


.1710?7^20 

.  1072^8610 

.050628216 

.023178447 

.011713617 

.006719022 

.0  04'J 

.003134f   • 

.002462f 

.002069167 

. 

.001671 • 

.00] * ■ 

I  -.98018 
.    1448363 

1-1  '651 


APPENDIX  F»  CONTINUED 


12H 


RUN  NUMBER  C-l 

rreak  frfquency 
damping  factor 
error  squared  = 


2.759 
.5643 
00016278 


FREQUENCY 


EXPERIMENTAL 


CALCULATED 


1.0000 

2.0000 

3.0000 

4.0000 

5.0  0  00 

6.0000 

7.0000 

8.0000 

9.0000 

lu.0000 

11.0000 

12.0000 

13.0000 

14.0000 

15.0000 

16.0000 


.140809650 
.132512380 
.089797529 
•033729233 
.010363286 
.005779744 
.003182099 
.002182841 
.001469301 
.001242401 
.001082352 
.000883831 
.000664090 
.000613566 
.000663485 
.000717838 


.135808600 
.139951850 
.081517003 
.032498836 
.013698786 
.006682771 
.003725341 
.002322848 
.001591469 
.001180310 
.000934858 
.000780990 
.000680553 
.000612726 
.000565576 
.000531972 


RUN  NUMRFR  C-2 
PREAK  FRFQUFNCY  =  2.491 
DAMPING  FACTOR  =  .6106 
ERROR  SQUARED  =  .00014516 


FREQUENCY 


EXPERIMENTAL 


CALCULATED 


1  .0000 
2.0000 
3.0  000 
A. 0000 
5.0  000 
6.0  000 
7.0000 

a. oooo 

v. 0000 

I  .0  000 

I I  .0  000 
12.0  000 

i  -i.nooo 
1  4 . 

1  '■  .0  000 
1 6.0000 


.162035830 
.13331 1480 
.073^1 5408 
.021932241 
.007773081 
.004792393 
.0024-13823 
.0015  72545 
.00  1 1  4  12V7 
•00089697< 
,000659970 
•00049839] 
.0004."  !   ■ 
,0004341 
,000448^ 
,00  04  7  077  7 


.  1  «SQS74390 
.1^8545480 
.06383m  ] 
,02  1980520 

•  0101'  ■  • 

.Oi  ;.'7 
.00;'  771417 

,001  M 

,001 148 

.  00831511 

,0006401 

•OOOS20738 

.0 

,01 

•  01 

. 


APP       f-  ,  CCNTIN 


:tcr 

-.naorn  =  .0000^^6 


FPFOUF"   ' 


FXPF PI"FNT ■ 


CALCULATED 


1  .  0 
2.0000 

.  noo 

•  .   00 
.  '000 

•  .   00 
7.0000 
8. 
9,0000 

l l .oooo 

00 
.     ~00 
14. 

1  '   . 
16.00CC 


.  18590 
.10  191 
.05010 
.01681 
.00658 
39C 
.0023*3 
.oniflf, 

,001 60 

,nop8r> 

.O^OAQ 

•00078 

.00074 
.00069 
.00062 
.00051 


2900 
1160 

6079 

9440 
6521 
3263 

9158 
41  ^9 
8^90 

8422 
1877 
11  12 
2646 
3682 


.  1852 
.1047 
.0430 
.0176 
.0083 

.0026 
.0017 
.0012 

.0009 
.0008 
.0006 
.0006 
.0005 
.0005 
.0005 


^^50 

42670 

80219 

64767 

73  A 

63098 

67067 

6049R 

I  22 

79000 
■86 
91244 
16982 
66258 
30672 
05119 


RUN  NU'-'PFR  P-l 
npp»<  FREQUENCY  =  ?.74c> 
FACTOR  =  ,7R18 
2Q  c^,,APpn  -  .00016406 


FREQUENCY 

EXPFRI  MENTA 

. 

.162367470 

.  "^000 

•106652920 

. 

•068475306 

4.0000 

. 

. 

.01 n^A0945 

. 

.00S919274 

. 

,001844248 

00 

, r     2442703 

9. 

•  0       2  8  6 

00 

.001265954 

i 1.0000 

.001092988 

. 

.  00807471 

1  . 

.000550864 

. 

•  0       7  59 

'  . 

497 

■  . 

860 

CALCULATED 


.  1601 
.  1136 
.0584 
.0269 

•  0130 
.00<S8 
.00^9 

•  0( 

.001  l 
.0008 
.0007 
.0005 

.0004 
.  )03 


3  2  770 
55940 
47364 
2^139 
1  ^794 
79875 
73884 
841  1  1 

rei 

05392 
80888 

93( 
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APPENDIX  F»  CONTINUED 


RUN  MUMBFR  E-l 
BREAK  FRFOUFNCY  =  2.397 
DAMPING  FACTOR  =  .8198 
ERROR  SQUARED  =  .00039070 


FREQUENCY  EXPERIMENTAL  CALCULATED 

1.0000  .176176300  .174257760 

2.0000  .093446826  .102137000 

3.0000  #061565092  .044405686 

4.0000  .021132733  .019011632 

5.0000  .005585339  .008993197 

6.0000  .004021728  .004757374 

7.0000  .002308439  .002782796 

8.0000  .001702735  .001777002 

9.0000  .001165651  .001225253 

10.0000  .000791363  .000903503 

11.0000  .000652180  .000706134 

12.0000  .000587770  .000579826 

13.0000  .000489721  .000496047 

14.0000  .000540654  .000438748 

15.0000  .000653968  .000398509 

16.0000  .000658421  .000369591 


RUN  NUMBFR  A-10-R 
BREAK  FREQUFNCY  =  2.552 
DAMPING  FACTOR  =  .7370 
ERROR  SQUARED  =  .00007221 


FREQUENCY  EXPERIMENTAL  CALCULATED 

1.0000  .166139080  .164948800 

2.0000  .11^465500  .117047270 

3.0000  .062388525  .055612244 

4.0000  .023223737  .023971784 

5.0000  .008370066  .011376313 

6.0000  .005547935  .006166417 

7.0000  .004704718  .003790209 

8.0000  .003548422  .002600341 

9.0000  .002375362  .001955926 

10.0000  .001812077  .001583736 

11.0000  .001545485  .001357092 

12.0000  .001^08950  .001212870 

13.0000  .001063812  .00111763* 

14.00CO  .000015109  .0010*7732 

15.0000  .000884'.  1'^  .00100728* 

16.0000  .000808838  .000O74701 


Apr        .    •  •  •  n 


-]-R 

■PING  ■  ,5432 

11540 


FREOU' 


EXPERIMENTAL 


CALCULATED 


1  . 

.  TOO 

. 
4. 

.  .000 
.  000 

1  .'000 

1 1 . 

12  . 

.   *00 
14. 


. 1 ? 90? 04 50 
.  •  /'QUARTO 
.0961  1  6r»61 
.  18710632 
,017657379 
,007542* 
,004507] 50 
14030 
.001964505 
.001554822 
.00] 170965 
.000906158 
,00081504] 
0903852 
,nop9^2895 

,0^0  8  44591 


.  1 . 

,  n  59  r,n]or> 

.090178471 

.017820799 

.01601 TP41 

.00782440? 

•0043966 

•002783738 

.001947676 

.001479887 

.001201642 

.00102770? 

.   .0914425 

.0O0818OSP 

.00078*048 

.00074711  2 


IBER  c-2-R 
A  K  FP        t  ■  2.607 
"PING  FACTOR  =  .6375 
29    SQUARED  =  .00016108 


FPf  ' 


FXPFRI'-'FNTAL 


r ALTULATPP 


1  . 
.   DO 

8. 

1  1  . 

■  . 


. l 562? 1 roo 

•125227890 

750 
•008713803 
.00  4  7  r.  8  8  08 
7  9  700 
.  i  ? 1064  06 
•001686 
. 

►02 


.  1  Sl?446^ 

.131601250 

•  06!  182 

. 

•01  587 

•  0058 
•00337 7861 

.0021 

.0      . 

.    1)0  9 1 5  < 
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ABSTRACT 

A  study  was  made  of  the  application  of  noise  analysis  to  the  determina- 
tion of  the  dynamic  behavior  of  a  copper-water  fluidized  bed  system.  The 
fluctuations  of  the  upper  fluid-bed  interface  were  recorded  as  a  function  of 
time  during  steady  state  operation  of  the  fluidized  bed.  The  input  to  the 
system  was  assumed  to  be  Gaussian  white  noise  and  the  system  was  assumed  to 
be  an  ergodic  random  process.  Under  these  conditions  the  significant  time 
constants  were  available  from  the  power  spectrum  of  the  system. 

The  theory  of  noise  analysis  was  summarized  with  particular  emphasis  on 
the  analysis  of  digital  data  spaced  at  equal  time  increments.  A  computer  pro- 
gram was  developed  to  calculate  the  autocorrelation  function  and  the  power 
spectrum  from  the  experimentally  observed  bed  height  data.  It  was  found  that 
the  power  spectrum  of  the  fluidized  bed  system  could  be  described  in  terms 
of  a  second  order  transfer  function.  Using  this  model,  the  break  frequency 
and  damping  factor  associated  with  the  transfer  function  were  determined  from 
the  experimental  power  spectrum  by  the  method  of  least  squares. 

The  break  frequencies  and  damping  factors  were  observed  to  vary  as  a 
function  of  particle  diameter  and  bed  expansion  ratio.  Plots  of  the  results 
indicated  that  the  break  frequency  decreased  as  a  function  of  bed  expansion 
ratio  while  the  damping  factor  increased  as  a  function  of  bed  expansion  ratio. 
No  such  trends  were  apparent  in  the  dependence  upon  the  particle  diameter, 
however,  statistical  analysis  showed  that  both  the  particle  diameter  and  the 
bed  expansion  ratio  had  a  significant  effect  on  the  break  frequency  and  damp- 
ing factor. 


Step  input  responses  were  recorded  and  the  experimental  results  were 
compared  with  those  calculated  from  the  second  order  transfer  function.  It 
was  found  that  the  results  compared  quite  well  for  both  positive  and  negative 
small  perturbations.  For  large  step  inputs  the  transfer  function  predicted 
more  overshoot  than  was  observed  experimentally.  The  results,  however,  in- 
dicated that  noise  analysis  should  not  be  overlooked  as  a  method  of  deter- 
mining the  significant  time  constants  of  a  fluidized  bed  system. 


